Ring exchange, the Bose metal, and bosonization in two dimensions 
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Motivated by the high-T c cuprates, we consider a model of bosonic Cooper pairs moving on a 
square lattice via ring exchange. We show that this model offers a natural middle ground between a 
conventional antiferromagnetic Mott insulator and the fully deconfined fractionalized phase which 
underlies the spin-charge separation scenario for high-T c superconductivity. We show that such 
ring models sustain a stable critical phase in two dimensions, the Bose metal. The Bose metal 
is a compressible state, with gapless but uncondensed boson and "vortex" excitations, power-law 
superconducting and charge-ordering correlations, and broad spectral functions. We characterize the 
Bose metal with the aid of an exact plaquette duality transformation, which motivates a universal 
low energy description of the Bose metal. This description is in terms of a pair of dual bosonic 
phase fields, and is a direct analog of the well-known one-dimensional bosonization approach. We 
verify the validity of the low energy description by numerical simulations of the ring model in its 
exact dual form. The relevance to the high-T c superconductors and a variety of extensions to other 
systems are discussed, including the bosonization of a two dimensional fermionic ring model. 



I. INTRODUCTION 

Despite intense experimental effort exploring the phase 
diagram of the cuprates, the nature of the pseudo-gap 
regime remains mysterious. In the very underdoped nor- 
mal state there are strong experimental|-Jjiints of local 
pairing and superconducting correlationsaau, despite the 
absence. of phase coherent superconductivity. Within this 
pictures, the pseudo-gap regime supports a pair field 
with non-zero amplitude but with strong phase fluctu- 
ations which disrupt the superconductivity. A theoreti- 
cal approach then requires disordering the superconduc- 
tivity by unbinding and proliferating vortices. Unfortu- 
nately, such an approach will likely face a most worrisome 
dilemma. Proliferation and condensation of sinsLs hc/2e 
vortices necessarily leads to a confined insulatoru, which 
should show sharp features in the electron spectral-func- 
tion in apparent conflict with ARPES experimental On 
the other hand, if pairs of hc/2e vortices condenseoQ, the 
pseudo-gap phase must necessarily supjaatt gapped un- 
paired vortex excitatipasQaa ("visons'^Ej'tj, which have 
yet to be observedOEl Is there an alternative possi- 
bility? In this paper we find and explore a truly novel 
quantum fluid phase of 2d bosons, which we refer to as 
a "Bose metal" , which answers this question in the affir- 
mative. The Bose metal phase supports gapless charge 
excitations, but is not superconducting - the Cooper pairs 
are not condensed. The hc/2e vortices are likewise gap- 
less and uncondensed. Being a stable quantum fluid with 
no broken spatial or internal symmetries, the Bose metal 
phase has many properties reminiscent of a 2d Fermi liq- 
uid - a Id locus of gapless excitations (a "Bose surface"), 
a finite compressibility, an almost T— linear specific heat 
with logarithmic corrections. 

Neutron scattering measurements in the undoped 
cuprates have revealed a zone boundary magnon disper- 
sion which can best be accounted for by presuming the 
presence of appreciable four-spin exchange processestJ. 



Recent theoretical work points to the importance of 
such ring exchange processes. .driving fractionalization 
and spin-charge separation&Ej'EEl. Motivated by these 
considerations, we focus on a class of model Hamiltonians 
describing a square lattice of 2d Bosons (Cooper pairs) 
with appreciable ring-exchange. It is hoped that such 
microscopic models will be appropriate in describing the 
charge sector in the under-doped cuprates. When the 
ring exchange processes are sufficiently strong, we find 
that the Bose metal phase is the stable quantum ground 
state. 

Lattice models of interacting bosons in two dimensions, 
such as the boson Hubbard model, have been st ud i ed 
quite exhaustively during the past several decadeslllllla, 
primarily as models for Josephson junction arrays and 
superconducting films but also in the context of quan- 
tum magnets with an easy plane U(l) symmetry and 
most recently in the context qf-trapped Bose condensates 
moving in an optical latticeO. Typically, the ground 
state phase diagram consists of a superconducting phase 
and one or more insulating statesEJ. At fractional bo- 
son densities commensurate with the lattice, say at half- 
filling, the insulating behavior is driven by a spontaneous 
breaking of translational symmetry. When accessed from 
the superconductor, such insulating states can be fruit- 
fully viewed as a condensation of elementary vorticeso. 
Very recent wack on a Kagome lattice boson model with 
ring exchangdia (which arises in the context of frus- 
trated quantum magnets) has revealed the existence of 
a more exotic fully gapped insulating state with no bro- 
ken symmetries whatsoever. In this state the charged 
excitations are-,"fractionalized" , carrying one-half of the 
bosons chargeO, and there are also gapped vortex excita- 
tions called visons. Here, we study the simpler system of 
bosons with ring exchange moving on a 2d square lattice. 
Our central finding is the existence of a critical gapless 
quantum phase of 2d bosons which is stable over a partic- 
ular but generic parameter range. This phase is in some 
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respects very similar to the familiar gapless Luttinger 
liquid phase of interacti:ag|-b|0S4)ns and Fermions moving 
in one spatial dimensionc3H-Eil. Indeed, to describe this 
phase we introduce a 2d generalization of bosonization. 
Specifically, we introduce a new duality mapping which 
transforms a square lattice model of 2d bosons (in a ro- 
tor representation) with ring exchange, into a 2d theory 
of "vortices" hopping on the sites of the 2d dual lattice. 
Then, based on the symmetries we construct a low en- 
ergy effective description, in terms the two (dual) phases 
of the boson and vortex creation operators. The Bose 
metal phase has a Gaussian fixed point description in 
terms of these fields. In addition, there are various non- 
quadratic interactions, which scale to zero in the Bose 
metal. Instabilities towards superconducting and insu- 
lating states are triggered when one or more of these 
interactions become relevant, closely analogous to 2fcp 
instabilities in a Id Luttinger liquid. 

Since this paper is quite long, the remainder of the 
introduction is devoted to a brief summary of the ring 
model Hamiltonian and the bosonized description of the 
Bose metal. The Hamiltonian we focus on describes 2d 
bosons in a rotor representation: 
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Here r labels sites of the 2d square lattice and cf> r and n r 
are canonically conjugate variables: 

[4>t, n r i] = i8 T>r >. (3) 



Representing the phase of the boson wavefunction, <j) r 
is taken to be 2-7T periodic <p r — <p r + 27r, so that the 
eigenvalues of the conjugate boson number operator n r 
are integers. The mean boson density is set by n, and 
we will primarily focus on the case with half-filling: n — 
i. The argument of the cosine in the second term is a 
lattice second derivative, involving the four sites around 
each elementary square plaquette. This term hops two 
bosons on opposite corners of a plaquette, clockwise (or 
counterclockwise) around the plaquette, and is an XY 
analog of the more familiar SU(2) invariant 4-site ring 



operators (which arises 



exchange term for spin one-ha 
in the context of solid 3 — HecBcB) . 

In addition to the conventional spatial, particle-hole, 
and total number conservation symmetries (see Sec. Ill), 
this ring Hamiltonian has an infinite set of other sym- 
metries. Specifically, the dynamics of the Hamiltonian 
conserves the number of bosons on each row and on each 
column of the 2d square lattice - a total of 2L associated 
symmetries for an L by L system. This is fewer than 
a gauge theory, which has an extensive number of local 
symmetries, but these symmetries will nevertheless play 
a crucial role in constraining the dynamics of the model 
and stabilizing a new phase. 



We will also consider adding a near neighbor Boson 
hopping term, 



«* = -* E 
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where the summation is taken over near-neighbor sites on 
the square lattice. This term breaks the 2L U(l) symme- 
tries corresponding to conserved boson numbers on each 
row and column, leaving only the globally conserved total 
boson number. 

The remarkable result of this paper is that, over a par- 
ticular but generic parameter range, ring models of this 
type sustain the novel Bose metal phase. The effective 
description of the Bose metal is in terms of low energy 
"coarse grained" variables cp ~ <f> and an additional field 
$ (see Sec. ill), which is related via n — n ~ Tr~ 1 A xy , d 



to the long-wavelength density, and further can be used 
to construct a "vortex creation operator" v ~ e 1 ® . To 
fix notation, we define the n r ,<j) T ,ip T operators on the 
sites of the original square lattice, with x, y coordinates 
taking half-integer values, and d r (with N r ,9 T operators 
to be defined later) on the dual square lattice with in- 
teger x, y coordinates. The ip, $ fields are governed by 
an approximately Gaussian (Euclidean) effective action 

S = J^drC, with C = Cq + C±, and r is imaginary time, 
(3 = \/k B T. The Gaussian part of the Lagrangian is 



Co = H [ip, 1?] + V* -dripr'Axy'dr, 
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where the sum is over sites r, r' = r + x/2 + y/2 on the 
dual and original lattices, respectively, and the effective 
Hamiltonian is 
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Here the momentum (k) integral is taken over the Bril- 
louin zone |&a;|,|fe y | < 7r, and (A xv (p)k denotes the 
Fourier transform of A xy (p r (and similarly for A xy -d). 
The functions /C(k),W(k) are non- vanishing finite peri- 
odic, and analytic. Their values along the k x and k y axes 
parameterize the Bose liquid, much like the effective mass 
and Fermi liquid parameters do a Fermi liquid. Experts 
will note a strong similarity to the bosonizj ed e ffective 
action for a one dimensional Luttinger liquid! 2 nri, which 
is explored in Sec. HI. Like in a Luttinger liquid, there 
are additional non-quadratic terms in the action. It is 
sufficient to keep only those of the form C\ — Ci v + Ci$, 
with 

£ « = EEC'" 008 ^^-^)], (7) 
r q,s 

with r on the original lattice, and 
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with r on the dual lattice. In both Eqs. Q — 2irn, 



r = (x,y), s = (m,n). Here, w™ q ' n = w^ m gives the 
amplitude to hop 2q vortices by the translation vector 
(771,71) (or (71,771)), and t™' n = t™' m gives the hopping 
amplitude for a q bosons along the same vector. In the 
Bose metal phase, all these terms are "irrelevant", i.e. 
give only perturbative corrections to physical properties, 
though these corrections can be significant. Also like in 
a Luttinger liquid, there is a non-trivial relation between 
microscopic quantities such as the energy and particle 
densities and the coarse grained variables of the low en- 
ergy theory. In particular, for the energy and particle 
densities, one finds 



l-q <-*xy u xy 
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where 8n = n — fi, e is the energy density, and c^ q are 
non-universal constants. 

The paper is organized as follows. In Section 2 we treat 
the ring exchange term within a "spin-wave" approxi- 
mation. This reduces the Hamiltonian to a quadratic 
form which can be readily diagonalized. Within the two- 
dimensional Brillouin zone there are gapless excitations 
along the lines k x — and k y = 0, which are associ- 
ated with an infinite set of conservation laws possessed 
by the ring exchange model. Dual "vortex" variables are 
then introduced via a new plaquette duality transforma- 
tion. This dual representation is well suited to numerical 
simulations since it is free of any "sign problems" . We 
implement a quantum Monte Carlo simulation, and show 
that the Bose metal phase is present over appreciable re- 
gions of the phase diagram of the ring-exchange rotor 
model. 

In Section III we construct the low energy effective 
model in terms of the dual vortex fields, and extract 
the universal properties of the Bose metal phase in Sec- 
tion IV. These properties are confirmed by the quantum 
Monte Carlo results. Section V is devoted to an analysis 
of the stability of the Bose metal phase. Remarkably, 
we find that for generic (incommensurate) boson densi- 
ties, there are regions of parameters where the Bose metal 
phase is stable towards all perturbations, even those that 
break the row/column symmetries. The situation is rem- 
iniscent of the Fermi liquid phase, whose "fixed point" 
description possesses an enormous set of "emergent sym- 
metries" - the number of Fermions on each patch of the 
Fermi surface are independently conserved. At commen- 
surate fillings we explore in some detail the instabilities 
of the Bose metal towards various insulating states. 

In Section 6 we consider a "hard-core" version of the 
boson Hamiltonian which allows us to obtain an exact 
zero energy wavefunction when the couplings are care- 
fully tuned. We then perturb away from the soluble 



point, and compute properties of the adjacent quan- 
tum phase which we thereby identify as the Bose metal. 
Finally, in Section 7 we conclude with a discussion of 
how various ring exchange processes emerge from a Z2 
gauge theory formulation of interacting electrons and 
then make some speculative remarks about the possible 
relevance of the Bose metal phase in the context of the 
under-doped cuprates. 



II. SPIN WAVES AND PLAQUETTE DUALITY 
A. Spin Waves, massless modes and symmetry 

Here we consider a spin wave approximation to the 
microscopic ring Hamiltonian which leads to a harmonic 
and soluble theory. The resulting energy spectrum van- 
ishes along two lines in momentum space. This remark- 
able feature is then shown to follow directly from the 
existence of the infinite set of conservation laws of the 
ring Hamiltonian. As we detail in later sections, such a 
harmonic description underlies an effective theory of the 
Bose-metal phase, in close analogy to bosonization of the 
Id Luttinger liquid. 

To this end, it is useful to employ a combination of path 
integral and Hamiltonian methods. A standard path in- 
tegral representation is constructed in the usual way us- 
ing <f> eigenstates. In the time continuum limit, the par- 
tition function for the pure ring Hamiltonian takes the 
form, 
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with the Lagrangian 



/Ji^j( d T<t>r) + ind T (f> r - K cos(A xy (f> r )}. (12) 



The resemblance of the above Lagra&gian with that for 
the standard Boson-Hubbard modelEZI (which has a sin- 
gle, rather than double, lattice derivative inside the co- 
sine) suggests that one might try expanding the cosine 
potential to quadratic order. Doing so gives a soluble 
Harmonic theory for the action, which can be readily di- 
agonalized as, 
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with k = (k x ,k y ) living in the first Brillouin zone, 
k x ,k y S [— 7r, tt]. Remarkably, the energy dispersion, 



E k = AVUK\ sin(k x /2) sin(fc y /2)|, 



(14) 



vanishes on both the k x and k y axis. The presence of 
these gapless excitations can be traced directly to the 
existence of the infinite set of symmetries of the ring 
Hamiltonian, which conserves the number of bosons on 
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each row and on each column of the 2d square lattice. 
Specifically, these symmetries imply an invariance of the 
energy (and action) under, 



(15) 



for arbitrary functions <fr x (x) and $ a (y) . This invariance 
dictates that in the Harmonic spin wave form the energy 
must vanish whenever k x = or k y = 0. 

Since the plaquette term involves a lattice second 
derivative rather than an ordinary lattice gradient, how- 
ever, it is clear that phase {<$) fluctuations will be large in 
the ring model, and the spin wave expansion is suspect. 
Indeed, even in the classical limit — > 0, one can read- 
ily see that "vortex" configurations in which <fi r winds by 
2n around some plaquette (e.g. for a vortex with center 
at x = y = 0, 4> xy = n/2e(-x)Q(y) + nQ(-x)e(-y) + 
3ir/2®(x)<d(~y), where &(z) is the Heavyside step func- 
tion) are Unite in energy rather than logarithmic as in 
an ordinary XY model. Further, "double vortex" con- 
figurations in which this winding is 4-7T can be smoothly 
deformed into zero energy conformations (e.g. 4> xy = 
7r0(— x) + 7r6(— y)). This suggests that for non-zero 
U there will be substantial "vorticity" in the low-lying 
states. 

To address the legitimacy of the spin wave expansion, 
it is necessary to account for the periodicity of the cosine 
potential. This is most readily accessed by transforming 
to a dual form, just as one transforms to a dual bosonized 
representation for a system of Id particles£-l Here we 
introduce a new 2d quantum duality transformationEH 
which is specifically tailored for the ring model. 



To interpret the new variables, consider now a vor- 
tex centered on some site (x, y) of the dual lattice (the 
"core" ) . Classically, for the four sites on plaquette of the 
original lattice surrounding this site, 



oi<t>. 
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where a, b = ±1, and N v denotes the number of vortices 
(vorticity) on this plaquette. Comparing to Eq. JiH}), 
one finds TV = N v (mod 2), so N can be interpreted as a 
vortex number operator, modulo two. Since 9 and N are 
canonically conjugate, the operators e^ tB perform canon- 
ical "translations" of N, and hence can be regarded as 
vortex creation and annihilation operators, respectively. 
Note that since 8 € nZ, e 2lB = 1, consistent with the 
ambiguity in N = N v under even integer shifts. Clearly 
periodicity of <f> in the original variables is encoded in the 
discreteness of 9 in the dual description. 

The plaquette Hamiltonian when re-expressed in the 
dual variables reads, 



U 



H a = -Kj2cos{nN r ) + —J2l& 



(20) 



We will for the time being neglect the tunneling term H t , 
which will be returned to later. Note the strong similar- 
ity between the dual and original forms. To bring this 
out more clearly, and for the numerical simulations to be 
considered shortly, it is useful to go to a path integral 
formulation. In particular, consider the partition func- 
tion 



B. Plaquette Duality 

We consider dual fields living on the sites of the dual 
lattice, denoted 9 r and N r , which are canonically conju- 
gate variables: 



[N r , 9*] = id r>1 



(16) 



Here we take 7r7V r to be 2ir periodic so that 9 t /tt has in- 
teger eigenvalues. The reason for this unusual choice of 
normalization should become apparent below. For nota- 
tional ease we are denoting the sites of both the original 
and dual square lattices as r. In analogy to Eq. (j^), 
it is convenient to introduce an operator defined on the 
plaquettes of the dual lattice (i.e. sites of the original 
lattice) as: 



A xy 6 r — r — 



J r+y 



(17) 



These two new dual fields are related to the original fields 
by the relations, 
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where r' = r + x/2 + y /2. One can check that with N r 
and 9 r as conjugate fields, the original variables satisfy 
the required commutation relations. 



Zn = Tr e 



-0H a 
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Expanding Zn in the usual Trotter fashion with a time- 
slice Ar = e — > + , using the discrete basis of eigenstates 
of 9 r , one finds (Appendix ^) 
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with the "Lagrangian" 
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where d T 9 r = [6 t (t + e) — # r (r)]/e. The e-dependence 
of the time derivative term in Eq. ( p3[ ) is familiar 
from the "time-continuum limit" relating e.g. the 
d + 1-dimensional classical and <i-dimensional quantum 
transverse-field Ising models. Note the strong similarity 
of Eq. (23) to Eq. (|l2|), which emphasizes their nearly 
self-dual nature. 

The formulation in Eqs. ( p2j]23] ) is quite convenient for 
numerical simulations. For the simulations, we define an 
integer- valued "height" field 9 t (t) through 



A IB 9 r (T) = 7m( s > + 7rA xy r (r) 



(24) 
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FIG. 1: The phase diagram of the time discretized dual model 
as a function of the two coupling constants, showing a com- 
pressible Bose metal phase and an incompressible phase which 
we identify as a (n, n) charge density wave. The compress- 
ibility « was obtained from the k — > extrapolation of the 
density susceptibility x nn (k, u) = 0) on 8 x 8 x 32 lattices 
as indicated, with some checks made on larger system sizes. 
As seen from the inset, the critical U / K stays finite in the 
time continuum limit eK — > indicating a phase transition 
at U/K « 2.4 ± 0.4 in the quantum model. 



where = — l) x+v ]/2 is the background staggered 
density, chosen such that 9 t (t) obeys periodic boundary 
conditions. In this language, we obtain a generalized 
(anisotropic) solid-on-solid model in 2 + 1 dimensions. 
The correlation functions for the ^-fields are easily re- 
expressed (and numerically evaluated) in terms of these 
height variables, using Monte Carlo methods. An indi- 
cation of the sort of results obtained in shown in Fig. |l| 
which presents the phase diagram of this model as a func- 
tion of U / K and e, based on an evaluation of the den- 
sity correlations. For the simulations, we worked on a 
L x xL x x L T lattice, with various system sizes indicated in 
the paper. We used a Metropolis algorithm with a single- 
site update, 9 t (t) — > 9 T (r) ± 1, which corresponds to a 
ring-exchange move for the boson density. We checked 
for equilibration of various quantities, and averaged the 
data over 10 4 -10 6 sweeps of the lattice depending on the 
correlation function and location in the phase diagram. 
Notice that in the time continuum limit (e — > 0), the 
simulations reveal two phases as the dimensionless ratio 
K/U is varied, separating an ordered charge density wave 
state at large U from the Bose metal phase when the ring 
exchange term is large. 

Before obtaining an effective low-energy description, 
it is convenient to rewrite the dual partition function, 
Eq. (B^) in terms of continuous 9 variables using the Pois- 



son summation formula, 

Z = f [d9 r (r)} £ e £>)(2- r M0 r (r)-^ r MI) e - £ £ T L n[e ] 

s r (r) 

(25) 

In Eq. (Eq), we have included a parameter S <C 1, which 
"softens" the discrete- 9 constraint (exactly imposed for 
5 = 0). Carrying out the sum over the integer- valued 
s f (t) variables, one finds 



-Sq— Si 



Z = Z J [t»r(r)]e" 
where Zq is a constant (divergent for 5 — > 0), 
So = 



Si = -^^V2 9 COS2g6» r (T) 
rr q—1 



(26) 

(27) 
(28) 



and the v<i q are 0(1) coefficients whose precise values are 
not important. It is convenient to drop the constant and 
introduce infinitesimal source fields h p , h £ , for the density 
and energy density, respectively, 



Z[h",h E ] = J [d9 r (r) 



— So— Si— S h 



(29) 



with 



S h = ~ }Z iK-ij) ( A xy°r - Trn) + K(T)(d T 9 r ) 2 ] , (30) 

r.r 

where r' = r + x/2 + y/2. 

III. EFFECTIVE MODEL 

We now make an educated but perhaps bold guess 
as to the nature (and existence!) of a low energy ef- 
fective description. In the spirit of the renormalization 
group and effective field theory, we imagine defining a 
temporally coarse-grained variable i? in which the high- 
frequency modes of 9 are averaged over 



i? r (r) = [9 r (T)] f - nnxy, 



(31) 



where the square brackets indicate an average over "fast" 
high-frequency modes, and we have for convenience re- 
moved the mean "curvature" in ■& by a non-fluctuating 
spatially dependent shift. Provided that we average only 
over high-frequency modes, we expect that the resulting 
effective description of $ will remain local both in space 
and (imaginary) time. Inspection of Eq. (|2fj| ) shows that 
the dual microscopic action for 9 can be written as a sum 
of a quadratic part (£ r r eCn) and non-quadratic correc- 
tions. We postulate that the low-energy effective action 
for $ is a sum of a renormalized quadratic term similar 
to Cn and small non-quadratic corrections. Although we 
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have not implemented a detailed renormalization group 
treatment, the latter statement is tantamount to declar- 
ing the system is controlled by a stable Gaussian fixed 
point. Obviously such a fixed point can have at best 
a large finite basin of attraction, and we will return to 
the problem of determining the region of stability of the 
Gaussian theory in Sec. M. 



• Four-fold rotations: 9 x>y — > 9 y - x , K+i/2, v +i/2 ~ ¥ 

— h p h £ — > h e 

"j/+l/2,-x-l/2' IL x,y ^ n y-x- 

• Particle-hole symmetry: This is an invariance only 
for integer or half-integer densities (2n £ Z). For 
such values, the symmetry operation is 8 xy — > 
2ifnxy — 6 xy , h p — > —h p , h £ r — ► h e T . 



A. Symmetries 

Formally, the above postulate of the effective field the- 
ory description can be formulated as the statement that 
the generating functional can be approximated for low- 
frequency source fields by Z[h p ,h £ ] as ZoZ[h p , h s ], where 
Zq is an irrelevant constant, and 



[dMr)}e- SoW - SlW - 



(32) 



where Sq [i?] is a renormalized quadratic form, <Si contains 
small renormalized non-quadratic perturbations, and Sh 
is linear in h p , h € (we drop higher order terms in the 
infinitesimal sources). In general, we expect all three 
parts of the action to take the most general form possi- 
ble consistent with locality, analyticity (since only high- 
frequency modes are averaged over), and the symmetries 
of the problem, which should be found from the micro- 
scopic partition function, Eq. (p9|). These are 



Row/Column symmetry: 9 x , y {j) — > 9 x , y ij) + 
ir[m x (x) + m y (y)], where m x (x) and rn y (y) are ar- 
bitrary integer- valued functions. This corresponds 
to conservation of the number of vortices (modulo 
2) on each row and column, a dual version of the 
conserved boson number on each row/column of the 
original lattice. 

Space and time translations: 6 t (t) — ► O t +-r.(t + tq), 
where R is an arbitrary lattice vector and tq is 
an arbitrary real number. Under this translation, 
the sources must also be translated, h p ^ e (r) — > 

Reflections across a row or column containing a site 
of the dual lattice (or bonds of the original lattice): 
'sx s'y, where s, sf = ±1. The sources 

hi 



® X ,y * ^ ^ 

transform as h p 



ss'h p , 

"'sx,s'y' 



"x,y 



sx,s'y 



• Reflections across a row or column containing 
bonds of the dual lattice (or sites of the orig- 
inal lattice): This is not independent, and can 
be obtained from a composition of a trans- 
lation and a site-reflection. But for com- 

ss '9l/2+s( X -l/2),l/2+s'(y-l/2)7 



pleteness, 
h 
h 



ss "l/2+s(a;-l/2),l/2+s'(a-l/2)' x,y 



l/2+s(x-l/2),l/2+s'(j/-l/2)- 

• Time-reversal: r (r) — s 

h P/£ (~T). 



,(-r), h p/£ (r) 



B. Gaussian action 

The form of S is dictated by these symmetries and 
the relation between the microscopic and coarse-grained 
fields, Eq. (31). Note that the shift by iffixy implies that 
•Q is not a scalar. Consider first Sq. By space and time 
translational symmetry, it is diagonal in momentum and 
frequency space, 



S ° = lJ k J M(k,^Mk,c)| 2 , 



(33) 



where / k = /d 2 k/(27r) 2 , J u = f^ckj/^), and the 
momentum integral is taken over the Brillouin zone 
\k x \, \k y \ < 7r. By row/column translational symmetry, 



M (k x — 0, k y , U) = 0) = for all k 



(34) 



and similarly for k x «-> k y . The latter condition implies 
the existence of gapless excitations along the k x and k y 
axes in momentum space. These are the only such gapless 
states mandated by symmetry, and we expect the low 
energy physics to be dominated therefore by momenta 
near these axes and by low frequencies. Quite generally, 
the kernel M. can be expanded near the zero modes at 
k x = u) = 0, and at lowest order takes the form, 



M(k,L)) = A{ky)u) 2 + B(ky)kl 



(35) 



where the expansion functions A(k y ) and B(k y ) are even 
and 27r-periodic in k y . A similar expansion is of course 
possible around k y = it) = 0, with the identical expansion 
coefficients due to rotational symmetry. Analyticity at 
k = further implies B(0) — 0. A convenient represen- 
tation of M. which satisfies these requirements is 



M(k,u) 



El 



^ 2 /C(k) 



with mode energy, 



4yX(k)Z7(k)| sm{k x /2) sin(fc y /2)|. 



(36) 



(37) 



Here, W(k) and /C(k) are positive 27r periodic functions 
which characterize the Bose metal phase, but it is only 
their behavior for k, 

universal low-energy properties of the theory. 



<C 1 or k y <C 1 that determine the 
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C. Interaction terms 

Consider next the interaction terms, Si — JdrCi. Lo- 
cality requires that they couple combinations of i9 r fields 
only with nearby points r. Hence we consider a succes- 
sive sequence of terms, C\ = ^ m £i|mi coupling d fields 
at a total of to distinct points. We expect that the C±. m 
becomes (exponentially) increasingly small with increas- 
ing to. First consider the single-site terms, which are 
highly constrained, particularly by translational invari- 
ance, under which d xy + irnxy transforms as a scalar. 
Generally, 



a ;1 = -E 



£ v 2q cos(2 g i? + qQxy) + /C" 1 {d T V) 2q 



xy Lg=i 



q=2 



(38) 



where Q = 2irh, and has the physical meaning of the 
smallest reciprocal lattice vector of a one dimensional lat- 
tice with density n (like 2kp for a charge density wave). 
The second set of terms involve more time-derivatives 
than the analogous quadratic interactions in Sq, and 
hence are negligible at low frequency, so we will take 
ICq 1 = in the following. The viq terms will play an 
important role in Sec. [V]. Note that, except when the 
density takes special commensurate values (such as the 
interesting case n = 1/2), they are strongly oscillatory in 
space. 

Next consider two-site terms. Similar to the K~ l 
terms in Eq. ( |38| ) , a variety of spatial (lattice) and time 
derivative invariants are possible, but are negligible rela- 
tive to 5o, so we do not include them here. More inter- 
esting are sinusoidal terms, which must take the form 



A;2 = - ^lw 2q2q , mn cos[2q^ xy ~2q''& x+mty+n + Q{qxy-q'{x + rn){y + n))} 

xy;m,n qq' ^ 

+ Q{qxy + q'(x + m)(y + n))\ 
I 



(39) 



where 
„..± 



rotational 
„..± 



2q,2q' ,nm 



we may also impose w^ q , 



and reflection invariance require 
and without loss of generality 
Like the v 2 Q 



2q,2q' ,mn> 



2q,nm 



W 



2q,2q' ,mn ' 



2q 

terms in Eq. (|38|), most of the operators in Eq. (|3£j) are 
highly oscillatory for generic n. An important exception 
arises in w^ q 2q mn for rational densities n = z' / z, where 
z, z' are integers. Then these terms are non-oscillatory 
if z'mq and z'nq are z times integers. We consider in 
particular the case z' = 1, for which to, q and n, q can be 
chosen as all possible factors of z£ into two integers, with 
arbitrary integer I. Keeping only these terms, one has 



£i; 2 = -*y] w 2 n q n cos[2q(& xy - 



x-\-m,y-\-n 



)],(40) 



xy m-q—n-q—zl 



'•2q 



2q 



by rotational invariance. Since e 



it) 



acts as a vortex creation operator, the terms in Eq. (|40|) 
can be thought of as hopping 2q vortices a distance s = 
(to, n) on the 2d dual lattice. 

Finally, consider Sh — JdrCh- As for L\, we expect 
that the largest contributions will occur in terms involv- 
ing "i? r fields at a small number of distinct points r. For 
simplicity, we keep only the most local of these. Employ- 
ing the symmetries above, one finds 



^ h _ X! I K+l/2,y+l 



1/2 



<--Q <-*xy u xy 



£<5. 

q=l 



COs{2qd xy + qQxy)} 1,(41) 



where c^/ q are constants determined by the high-energy 
physics. Here we have kept only the single- plaquette 
terms in the charge density, and the single-site terms in 
the energy density. Note that the operators multiplying 
h p involve sines rather than cosines, which is a conse- 
quence of particle-hole and reflection symmetries. The 
form in Eq. ( [i"l|) can be derived perturbatively in 
using an explicit one-step coarse-graining proceduretil. 
The quantitative reliability of this constructive procedure 
is, however, limited to small bare V2 qi so for the prob- 
lem of interest the constants c^ £ should be viewed as 
phenomenological parameters. Eq. (^Tj) should be inter- 
preted analogously to the non-trivial bosonization rules 
for e.g. the charge density and other operators in one 
dimension. Very generally, a low energy effective action 
must be accompanied by rules specifying the operator 
content of various physical fields. In analogy to the no- 
tations used there, one may rewrite Eq. (|ll]) as a scal- 
ing equality between the microscopic and coarse-grained 
fields, as in Eq. (Q). 



D. Self-dual form - Bosonization analogy 



X] c 2q^xy sin{2qd xy + qQxy) 

9=1 



Having established the form of the effective action for 
the -d variables, we are also interested in being able to 
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compute quantities involving the boson phase ip. Since 
the effective theory is perturbative in <Si, it is straight- 
forward to re-introduce (p using Gaussian integration. 
In particular, the quadratic action 1S0 in Eq. (j33|) with 
the kernel in Eq. d36| ) can be transformed into the form 
shown in Eq. (|5j) in the introduction using a Hubbard- 
Stratonovich transformation: 

J[d V ]e-^ L j: r i^ r ,A xy9 ^ (42) 

where r' = r + x/2 + y/2 as usual, and (A xy ip)k denotes 
the Fourier transform of A xy tp r . After an integration by 
parts in the last term on the right hand side above, the 
full Gaussian part of the action takes a particularly trans- 
parent form, So — J Co, with Lagrangian Co defined in 
Eqs. di). 

The partition function is now represented as a path in- 
tegral over both sets of "low energy" fields, ip r and i? r , in 
an appealing self-dual form. The Gaussian theory above 
gives a fixed point description of the Bose metal phase. 
When augmented by the operator content of the fields, as 
in Eq. (0), together with the irrelevant non- linear inter- 
action terms in Eqs. ( (38|]3S| ), it gives a complete descrip- 
tion of the universal properties of the Bose metal phase, 
as detailed in the next Section. 

It will sometimes be convenient to integrate out the 
dual field, leaving a Gaussian action just in terms of 
the phase of the boson wavefunction: 

S ° = ^LI M <pOt,u)\<P(&u)\ a , (43) 
with Kernel, 

and with as given in Eq. ( ^7|) . Notice that this form 
is identical to the Gaussian fixed point theory in terms 
of d in Eq. (|33|), except that U(k) has replaced 7r 2 /C(k) 
in the denominator of Eq. ( ]44| ) . 

The above self-dual representation makes particularly 
apparent the close analogy betweenj— this theory and 
Bosonization theory in one-dimensiorjEj. In particular, 
underlying both theories is a pair of dual scaler fields, one 
the phase of the wavefunction and the other related to 
the density. The Luttingcr liquid fixed point is a Gaus- 
sian theory in the two fields, and has a form which is 
nearly identical to the Bose metal action above: 

Cid = H ld + - [ d T 4>d x 9, (45) 

7T Jx 

with Gaussian Hamiltonian: 

Hld = L [ Y {d ^ )2 + ^ {d ' J)2] - (46) 



Notice that the long-wavelength particle density in 
Bosonization theory, n ~ d x 6/-K, has simply been re- 
placed by a (lattice) second derivative in our 2 + 1- 
dimensional theory: n ~ A xy d/-K. The commutation 
relations between the two dual fields has also been modi- 
fied, as is apparent by inspecting the Berry's phase term 
involving both fields in the above Lagrangians. The ex- 
pressions in Eq. (^|) relating the bare boson density and 
energy density to the low-energy fields are a generaliza- 
tion of the more familiar Bosonization expressions, where 
for example the Boson density near 2kp is proportional 
to cos(26»). 

A key strength of Bosonization in Id is that it allows 
one to study the instabilities of the Luttinger liquid to- 
wards various types of ordered phases, such as a charge 
density wave state. Similarly, the above Gaussian rep- 
resentation of the Bose metal fixed point is particularly 
suitable for studying instabilities both towards insulating 
states with broken translational symmetry and towards 
a superfluid. But before undertaking this analysis, we 
explore the Bose metal phase in some detail. 

IV. THE BOSE METAL PHASE 

We now turn to the physical properties of the Bose 
metal phase, using the effective low energy theory devel- 
oped in the previous section. For the present, we will 
assume stability of the Bose metal, so that all physical 
quantities are accessible via perturbation theory (in Si) 
around the Gaussian theory for cp. The validity of this 
assumption is discussed in Sec. [v| 

A. Thermodynamic properties 

Consider first the thermodynamic properties of the 
Bose metal. The simplest is the compressibility k — 
dn/dfi. This is trivially given from the Gaussian the- 
ory, and one finds k = l/U(k = 0). Numerically, the 
compressibility may be obtained from an extrapolation 
of the density-density correlation function, 

X nn(K^n)=J2 J dr ^ (A xy 6 r (T) A xy 6 o (0)) e^ T -^ , 

from which k = Xn«(k — * 0,w = 0). The compressibility 
k determined in this way from the simulations were used 
to define the phase diagram in Fig. 

Also interesting is the specific heat. Since the Bose 
metal is essentially a free boson theory, this is determined 
completely from the density of states for the collective 
free boson modes 

p{E) = J S(E - Ek), (48) 
where the boson energy is 

£ k = tt(k)4|sin^-sin^|, (49) 
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and n(k) = v/W(k)/C(k). Like an ordinary (Fermi liquid) 
metal, the density of states gets a large finite contribution 
from the low energy modes near the "Bose surface" along 
the coordinate axes in momentum space. Unlike a Fermi 
liquid, however, there is a weak logarithmic divergence 
associated with the crossing point at k x = k y = 0. In 
particular, 



p(E) 



[ln(l/e)+C + O(e)], 



where Sl = f2(0, 0), e = E/$l and 



(50) 



with 



K(k) 
2E V 



(1 



(57) 

Since Ek vanishes linearly both as k x — > and as k y — ► 0, 
inspection of Eq. ( |57|) shows that F^(t) has a double log- 
arithmic divergence as r — > oo. Indeed, a careful calcu- 
lation shows that 



|lnr!or| 2 + Ffln(Oor) 



(58) 



Co = 41n2 



dk 



o sinfc/2 \Cl(k,0) 



n 



- 1 



(51) 



Note that the constant term in the density of states de- 
pends upon the full form of the dispersion f2(fc,0) = 
fi(0, k) all along the Bose surface, while the logarithmic 
term depends only upon the behavior at k = 0. The spe- 
cific heat is then 



c T - 



dx 



x 2 p{Tx) 
sinh cc/2 ' 

which to the same accuracy gives 

T 



C'i, 



47T 2 f7 



C /o + Ii 



(52) 



(53) 



and I n = f™dxx 2 (ln{l/x)) n /smh 2 (x/2), or J = 4tt 2 /3, 
7i = 2(3 - 2 7£; )7r 2 /3 + 8C'(2) w 4.643, where £(z) is 
the Riemann zeta function, and ~ 0.5772 is Euler's 
constant. 



for fl T 3> 1, with 



F[ = 




(59) 



[21n7T + 7 jE ], 



(60) 



where U a = U(0, 0), fC = /C(0, 0). The In 2 n a r behavior 
of F^ (t) at large r implies that G$ (r) decays faster than 
any power law. This translates into a similar singular 
behavior for the tunneling density of states. Writing the 
boson Green's function in a spectral representation, 



/•OO 

G (O,r)= / dE Ptun (E)e- E ^, 
Jo 



(61) 



one finds using a simple saddle-point analysis that the 
above behavior at large r requires 



B. Boson correlation functions 

It is interesting to contrast the large density of states 
for collective excitations with the tunneling density of 
states for the original bosons, which we will see is strongly 
suppressed. In particular, consider 



G (r,r) 



/ e <*r(T)g- 



#o(0)\ 



(54) 



By symmetry, since the boson number is conserved on 
each row and column, G<f,(r, r) can be non-zero only for 
r = 0. To further determine the behavior of G</> in the 
Bose metal, we must relate the microscopic Bose cre- 
ation/annihilation operators to the low energy modes. 
By symmetry, 



e Wr [j + ,4cos(2tf r + 2nnxy) H ] 



(55) 



We expect the low energy properties to be dominated by 
the first term, i.e. simply replacing cj> — > ip. Since tp is a 
Gaussian variable, 



G^(r,r) 



(56) 



Ptun(E) ~ exp 



2 E 



l n ^ (/3mm ^ + A) 

E yh E ' 



where a 



(62) 

13 = F*, \ = F% — 1 — F$ + \b.F$. 
Thus, although the Bose metal possesses a large set 
of gapless modes, the density of states for adding a bo- 
son into the system vanishes at low energy. Again, this 
is analogous to a Luttingcr liquid in lDc-3; the conserva- 
tion of particles per row or column means that the added 
particle affects all the particles in a particular row or col- 
umn, leading to a suppressed amplitude for such tunnel- 
ing events. Remarkably, the tunneling density of states 
actually vanishes faster than in a Luttinger liquid, indeed 
faster than any power law in energy. This behavior can 
be roughly understood as arising from two orthogonal- 
ity catastrophes occurring simultaneously in the row and 
column in which the boson is added or removed. 

It is also instructive to consider the boson four-point 
correlation function, for simplicity at equal times. Due 
to the row/column symmetries this is non- vanishing only 
when the four points sit at the corners of a rectangle: 



(63) 
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Upon replacing <f> with the low energy field ip, this can be 
readily evaluated using the Gaussian Bosc metal action. 
For any fixed y, one finds power-law behavior in x (and 



vice versa since G^\y,x) = G£'(x,y)) 



with 



Gf{x,y) 



v(y) 



\v{y) 



For |a;| > 1. 



TV 2 I 

n JO 



dk 



IU(0,k) sin 2 (fct//2) 
/C(0,fc) sin(A;/2) 



(64) 



(65) 



Note that Eq. (H) requires only \x\ ^> 1 and places no 
restriction on y. Hence it obtains also when both argu- 
ments are large, and thus 



U/K=1.0 m 
eK=1.0 j 
L =L 




In t 



"^T\/ %%(lnx)(\ny)-C(\nx+lny) 



(66) 




fC(0,k) 2sinfc/2 \ K. k ' { ' 



C. Vortex correlation functions 

The exact plaquette duality makes it possible to define 
a number of characteristic "vortex" correlators in terms 
of the field e 10 . As for the boson correlator above, the 
behavior in the Bose metal depends upon the expression 
for the vortex operator in the low energy variables. By 
symmetry, we expect the leading terms (involving the 
smallest exponentials of ■&) to be 



Re [A 



(68) 



where in general A is complex. Note that, since 9/n is an 
integer, the microscopic vortex field satisfies e l9 — e~ l9 , 
but this is not true for In the special case of n = 1/2, 
particle- hole symmetry (i? — > — 1?) further implies A is 
real, whence 



,1/2 2Acos(i9 r - -xy). 



(69) 



With this relation in hand, let us consider first the 
vortex two-point function, 



G e (r,r) 



e ie r (.T) e -ie (.o)y 



(70) 



As for the boson correlator G$, the vortex two-point 
functions vanishes unless r — due to the dual row and 
column symmetries. Evaluating it using Eq. (pq) and the 



FIG. 2: Vortex correlation function Ge(0,r) = 
(e ier(T) e~ ie ° (0) ) in the Bose metal phase. The cor- 
relation function decays faster than a power law 

~ exp I m Qo T \ 2 ^j at long times (see text). In- 

set shows this decay in a finite size scaling plot (L = 4 to 
L — 32) of lnGe(0,r) versus lnr for r = L T /2, from which 
we extract UP- ps 1, nearly its bare value for these parameters. 



Gaussian action, one finds similar results to the boson 
tunneling density of states, 

G e (0,r)~exp(-^|lnfi r| 2 -F 1 £ '|lnr!oT|l, (71) 



where 



F. 



Ft = 




(72) 



/Co 1 
Uq k 



(73) 



In fact, these coefficients are obtained directly from 
Eq. (59) using the duality transformation /C(k) 



U (k) pit*. In Fig. || we show numerical results for the two- 
point vortex correlation function obtained from the quan- 
tum Monte Carlo simulation in the parameter regime cor- 
responding to the Bose metal phase. The downward cur- 
vature of the data is consistent with a decay more rapid 
than a power law, and as shown in the inset can be fit to 
the form in Eq. @ with %P- sw 1. 

Next consider the vortex four-point function, which 
due to the dual row/column symmetry is non- vanishing 
only when the four points sit at the corners of a rectangle: 

G$\x,y) =(e i6 ->{r) e ie,o{r) e ie oy ( T ) e ie xy {r)y (74) 
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This can be re-expressed in terms of the low-energy field, 
using Eq. (|68|), and then evaluated with the Gaussian 
Bose metal action. For any fixed y, one finds power-law 
behavior in x, 



G { e\x,y) 



cos(irnxy) 
\x\Vv(v) 



for Id > 1. 



with 



Vv (y) = J dk 



/AC(0,fc) sin 2 (fcy/2) 
U(0,k) sin(fc/2) 



(75) 



(76) 



Interestingly, one can see from Eq. ( |75|) that for the case 
n = 1/2, the four boson correlator vanishes exactly when- 
ever xy is odd as a consequence of particle-hole symme- 
try. This behavior, and the associated power-law cor- 
relations, are shown in Fig. ||. In the limit when both 
x, y — > oo, the vortex four-point correlator vanishes faster 
than any power law: 



(x, y) ~ cos(7ma;y)e 
as x, y — > oo, with 



%(lnx)(lny)-C v (Ux+lny) 



(77) 




With the exception of the cos(irnxy) prefactor in 
Eq. (|77j), all the results in this subsection can be obtain 
from those in the previous one by the duality transfor- 
mation ip <-y d and £/(k) <-* 7r 2 /C(k). 



D. Collective correlation functions 

Next, let us consider the correlations of "two-particle" 
operators such as the boson and energy densities. As 
above, we require the operator correspondences between 
the microscopic and effective variables. These were al- 
ready worked out in Sec. II], and summarized in Eqs. (^|- 
|Iot ) in the introduction. For simplicity, we focus on the 
most interesting case of ri = 1/2, for which Q — it, and 
furthermore keep only the lowest non-trivial harmonics 
with 9=1. In this limit, Eqs. (|j^|lO|) reduce to 



+4 



J- (_l)(*+a)(v+6)+o+& sin{2 $ x+a , y+b ) (79) 



a, 6=0,1 



c^ + c|(-l)^co S (2^). 



(80) 



Consider first the density-density correlation function, 



X„„(r,r) = {8n x+ i y+ i{T)8n± i(0)). 



(81) 




FIG. 3: The 4-point vortex correlation function Gg(x = 
l,y,r = 0), defined in the text, evaluated in the Bose metal 
phase. The correlation function vanishes for odd values of 
y due to particle-hole symmetry at n = 1/2. The expected 
power law decay of the envelope is shown in the inset as a 
finite size scaling plot, from which we extract the exponent 
T) v (l) ~ 1.9 ± 0.1, close to its bare value for U/K = 1. 



Substituting for 5n using Eq. (^9[) , one obtains three con- 
tributions to Xnn'- 



Xnn(r,r) 



(82) 

The cross term xZ is negligible. The first contribution, 
Xnn( r i T ) is just the correlator between "coarse-grained" 
densities cx A xy i} xy . This term is non-zero in the Gaus- 
sian theory, and gives a smooth function of r, with a 
power-law behavior at large arguments. For instance, at 
equal times and large |x|, \y\ S> 1, 



1 



Uq x 2 y 2 

More generally, Xnn nas a smooth Fourier transform, 



(83) 



Xnn (q,Wn) = 



El 



W(k) ul + E* 



(84) 



Perturbative corrections from <Si to Xnn do not modify 
this qualitative behavior. 

The remaining contribution to the density-density cor- 
relator comes from the sin 2$ terms in Eq. ( [f9| ) . Naively, 
this contribution is ultra- local (and hence uninteresting), 
i.e. vanishes unless |y| < 1, as a consequence of the 
fact that the discrete row/column symmetries arc pro- 
moted to continuous ones at the Gaussian level. One 
may interpret this as meaning that the "vorticity" on 
each row or column is conserved exactly in the Gaussian 
model. This conclusion, however, is incorrect once the 
non-quadratic corrections in S± are taken into account, 
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since expanding factors of V2(— l) xy cos 2§ can "supply" 
vorticity in units of 2 to a particular site. Hence, Xrm 
becomes non-trivial at second order in u 2 . Provided 
|a;|,|y| > 1, the appropriate second order perturbative 
expression is 



v 22 
Ann 



( r ~\ — (_-i\(x+c-a)(y+d-b)+a+b+ 
V ' / 8 l^abcd=0.1 \ L > 



c+d 



XC 



(4) 

x-\-c— a.y-\-d— b 



(85) 



where 



3 2i[tf 00 (0) + tf XH (T)-# :c o(Tl)-tfo l ,(T2)] 





(86) 

where the expectation value indicates an average cal- 
culated in the Gaussian theory. At large distances, 
\x\, \y\ > 1, one expects C^l c _ ay+d _ b {T) » C^(r), in- 
dependent of a, b, c, d. In this limit, only the prefactor in 
Eq. (|8^) depends upon a,b,c, d, and the sums can be car- 
ried out explicitly. To do so, it is convenient to employ 
the identity 

{-iyy = i[i + (-1)- + (-i)y - (-iy+y], (87) 



valid for integer x,y. Applying this identity to Eq. (85), 
only the last term survives the sum, and gives 



(88) 



for \y\ 3> 1. This indicates the presence of (ir,ir) 
correlations in the boson density. More generally, if only 
\x\ > 1 but \y\ > 1 but still of O(l), one finds 
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(r,r) £ il-(-l) v+d - b )C^ +d _ b (rl 



bd=0A 



. (89) 

so at any fixed y the density-density correlator oscillates 
at wavevector i as a function of x (and vice- versa by 
rotational symmetry). To establish the range of these 



correlations, we must consider C x . v (t) in some detail. 

(90) 



Using properties of Gaussian fields, we have 
Cl 4 y( T ) = Jdr 1 dr 2 exp [-c x ^(r, n, t 2 )] 



where 



(.r,Ti,T 2 ) 



(91) 



2 ( ($xy (r ) + oo (0) - $ x0 (ri ) - tf 0y (r 2 ) ) ' 



The calculation and analysis of c xy is somewhat in- 
volved, and is described in detail in Appendix [b|. There 
we derive a useful approximation for c xy which captures 
the behavior in all the relevant limits (\x\ 3> l,f2o|i" — 
Ti|,ftoN, and \y\ < M): 



2r lv {y)ln{x' + n 2 y T 2 ) 
+ yUoV [ \y\ 



^oN 



, (92) 



where fl y ~ f2o (but in general depends weakly upon the 
full form of K.(0,ky),U(0,ky), and y) and the crossover 
function !F(X) satisfies 



X< 1 
\nX X > 1 

|||0|), we find 

2 



(93) 



2 2\~ 2 Vv{y) 



Using Eqs. (B: 

£?W(r) w>i {£) "e 2 ^ (x' + ny) 

(94) 

and of course the same behavior with x <-> y. 

This power-law behavior of C^y , and hence x„„(r,r) 
translates into singularities in the static structure factor, 

S nn (k) = J2 e ~ Z{k * X+kvV) Xnn(r,T = 0), (95) 

a: J/ 

for wavevectors near 7r = (7r, it). In particular, let k = 
TV + q, and consider the limit <C 1 with fixed. In 
this limit the singular behavior of the structure factor is 
dominated by large \x\ but y of 0(1). Hence we may 
apply Eq. (|8^) to write 

Snniiz + q) iJki (96) 
W f J2 cos ft/fC 1 + cos %) E e <toX C$ (0). 

odd y x 

From the power-law behavior in Eq. (Q), one can readily 
see that the Fourier transform in x leads to singularities 
for small q x . Indeed, if any 7] v (y) < 1/4 for any odd 
y , Xnn diverges as q x — > at fixed g^, while for 1/4 < 
?7„ < 3/4, the structure factor remains finite but has a 
divergent second derivative at q x = 0: 

S„„(7r + q) ~ ^° m (q) + ^^( (Zy )sgn(a y )| fe |- Q «, (97) 

odd y 

where S„ n (q) is a smooth function, a y — 1 —4r) v (y), and 
^y(<7jf) ^ -u 2 COS( ?y2/(l + cosq y ) is a positive amplitude 
peaked at q y = 0. Hence the behavior for small q x is 
dominated by the minimum (over odd y) value of r] v (y) 
(maximum of a y ). 

The zero frequency susceptibility, XnnQ*) = 
X„ n (k, o;„ = 0) has a similar but stronger divergence 
due to the extra time integration: 



X«n(7T+q) ~ X°n (q) + /J A y <A/)sgn(%) |fe 



odd y 



(98) 



where X^„(q) is another smooth function, AP{q y ) oc 
v4^(<7 y ), and a y = 2 — An v (y) = l + a y signals the stronger 
divergence. The difference in exponents implies that if 
Snn(ir + q) has a divergent slope at q x — 0, Xnn{^ + q) 
itself divergences there. Conversely, if Xnn (t + q) has a 
slope divergence at q x — (occurs for 1/2 < n v < 3/4), 
the static structure factor does not. Numerical results 
for the density susceptibility in the Bose metal phase are 
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FIG. 4: The density susceptibility x«n(k, ui = 0) in the Bose 
metal phase, close to the phase boundary, along (0, 0) — > 
(27r, 27r). The cusp at (tt, n) is due to subdominant power law 
correlations in the Bose metal phase which may be viewed 
as arising from corrections to the density operator in the 
Gaussian theory as discussed in the paper. Inset shows the 
grayscale plot of the susceptibility over the entire Brillouin 
zone, centered on (n, n) (dark regions indicate large suscepti- 
bility, and the k x = and k y — lines are omitted from the 
grayscale plot for clarity). The singular cusp is visible along 
the indicated (0, 0) — > (2n, 2n) direction. 




(units of 



FIG. 5: The energy susceptibility X£e(k,c<J = 0) in the Bose 
metal phase, close to the phase boundary, along (0, 0) — > 
(0,2n). The cusp at (0,7r) is due to subdominant power law 
correlations in the Bose metal phase which may be viewed 
as arising from corrections to the energy-density operator in 
the Gaussian theory as discussed in the paper. Inset shows 
the grayscale plot of the susceptibility over the entire Bril- 
louin zone, centered on (n, tt) (dark regions indicating large 
susceptibility). The singular features are most clearly visible 
as peaks (dark spots) at the (0, it) and (n, 0) points and as a 
dip (white spot) near (7r,7r). 



shown in Fig. and reveal a singular cusp- like behavior 
at wavevector tt. This form is consistent with that pre- 
dicted by Eq. 
0. 

Now consider the energy-energy correlation function 



with a maximum value — 1 < a™ ax < 



X e£ (r,r) = (e xy (r)e xy (0)) c , 



(99) 



where the c subscript indicates the cumulant expecta- 
tion value. As for the density-density correlator, we can 
employ Eq. (|8Cj) to expand \ee as 



Xe£ (r,r) ~ (c^^(r,r) + ( C ^^(r,r)+^x^(r,r). 

(100) 

As for the density case, Xes nas a smooth Fourier trans- 
form, and Xee i s negligible. We focus on x^, which (for 
\x\, \y\ > 0) to second order in v 2 is 



(101) 



Using Eq. (g7|9J), one straightforwardly sees that there 
are singularities in Xee (k) as k x (k y ) — * and k x (k y ) — > n. 
In particular, 

S ee (k) ]kx ^S° es (k)+J2Al(k y )sgn(a y )\k x \- a y, (102) 

even y 

X £e (k)i^K< 1 xL(k)+^^(fc 2/ )sgn(5,)|fc :c r 5 «, (103) 

even y 



and 



S se (ir + q x , k y ) i te T«i (104) 
S° ee (n + q x , ky) + J2 A l ( k v) sgnK)|fe I , 

odd y 

Xe £ (7r + q x , k y ) i te T«i (105) 
X° £ (tt + q x , k v ) + Y,A e y {k y ) sgn(a v )\q x \- & « , 

odd y 

where A e y (k y ) oc A y (k y ) oc cos k y y. Similar formulae for 
k y w 0,7r are obtained by rotation. Note that because 
A y (k y ) tx cos k y y is negative for k y « tt and y odd, the 
energy-energy correlator has a singular dip near k = tv 
and singular peaks near k = (0, 0), (tt, 0), and (0, tt). Nu- 
merical results for the energy susceptibility in the Bose 
metal phase are shown in Fig. [|. Notice the singular 
cusp-like peak s at wavevectors (0,7r) and (tt, 0), as pre- 
dicted in Eq. ( |og ). 



E. Electrical Conductivity 

We finally consider the electrical conductivity in the 
Bose metal. To obtain an expression for the current 
operator requires coupling in a vector potential. The 
usual minimal coupling prescription for a lattice model 
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of bosons would have one replace a discrete lattice deriva- 
tive as, 



(106) 



with fi = x, y, and define the current operator by dif- 
ferentiating with respect to the vector potential. For the 
boson ring model this prescription is ambiguous due to 
the second derivative form, and a family of gauge inequiv- 
alent forms are possible: 



A 



xy^rr 



A 



xy^rr 



+ aA x Av + pA v A*, 



(107) 



where A x f r = / r +i — fr denotes a discrete derivative. 
Gauge invariance requires a + (3 — 1, whence this can be 
re-expressed as, 



A, 



A xy <j) r + A x A v r + (a- 1)$, 



(108) 



where $ = A X A V — A y A x is the gauge invariant flux 
through the plaquette. If one derives the boson ring 
model by starting with a model of electrons reformu- 
lated in terms of a Z 2 gauge field coupled to spinons 
and chargons as detailed in Appendix [S| one arrives 
at an appealing symmetric form for the ring term with 
a = P= 1/2. 

Henceforth we focus on the zero wavevector conductiv- 
ity, cr(u>). In this case, the above ambiguity is irrelevant, 
since a spatially uniform vector potential does not enter 
for any value of a. But the vector potential will still of 
course enter into the boson hopping term in Eq. (Q), and 
upon differentiation generates the usual current operator, 



K = isin( 



J r+fi 



(109) 



When we coarse grain the theory, we should write down 
the current operator in terms of the slow field, ip: 



1^ = C[t sin((y9 



r+fi 



fr) + ■■ 



(110) 



with c/ a dimensionless constant. The other contribu- 
tions will include terms which hop a single boson several 
lattice spacings and terms which hop several bosons - 
generally all local terms allowed by the symmetries. Such 
terms will generically be subdominant at low frequencies, 
and so we retain only the leading contribution. 

With the current operator in hand it is straightforward 
to obtain the conductivity from the usual Kubo expres- 



sion: 



with 



<7i(w) = Rea(iui n — > u>), 



(111) 



v(iu n ) = — f drJ2e- iUnT W(r)7g(0)). (112) 

We have dropped the diamagnetic contribution since it 
will not contribute to the real part of the conductivity. 
Evaluating the correlator using the Gaussian Bose metal 
action gives at zero temperature, 

(^(r)/ a (0))o ~ t 2 S yfi (x 2 + t 2 )- a , (113) 



with scaling dimension A = ry(l)/2 where rj{y) is given 
explicitly in terms of the Bose liquid parameters in 
Eq. 65 , Performing the time integration and spatial sum- 



mation gives a power law singular contribution, 

W„(7(iWn) = -At 2 {-lf^^- 2 + X reg (iu n ), (114) 

where A is a positive constant and S reg is analytic in 
it's argument and thus does not contribute to the real 
part of the conductivity. Analytic continuation to real 
frequencies gives a singular contribution to the complex 
conductivity of the form, 



a(u) = At^- 1 ) + w 3 (u), (115) 

with real 02 (w). We note that, causality places strong 
constraints on the phase anglec3 from the singular con- 
tribution, consistent with the above form. Finally, taking 
the real part gives the optical conductivity, 



ai(oj) = At 2 \ sin(7rA)||o; 



2A-3 



(116) 



Notice that the amplitude of this contribution vanishes 
for integer A. For these special cases, the higher order 
contributions to the current operator should be kept, and 
will contribute a similar form but generally with larger 
scaling dimensions. 

As we discuss in the next section, stability of the Bose 
metal phase to boson hopping requires that A > 2, im- 
plying an optical conductivity vanishing rapidly at low 
frequencies, <7i(w) ~ io a with a > 1. For A < 2, the 
Bose metal will be unstable to superconductivity, but for 
small boson hopping amplitude the transition tempera- 
ture would be low. In that case, the optical conductivity 
above T c might still be well described by the above power 
law form. 

One may also directly calculate the non-linear dc con- 
ductivity, <j(E y ,T) — d(Iy)/dE y as a function of E y and 
T. This is mathematically complicated, but formally 
quite similar to the perturbative calculations of tunnel- 
ing conductance between parallel Luttinger liquids, as 
carried out e.g. in Ref. We forgo this calculation 
here for the sake of brevity, quoting only the resulting 
scaling form 



cr(E y ,T) = At 2 T 2A - 3 G(E y /T) 



(117) 



where -> 1 as £ -> 0, and T{£) - £ 2A ~ 3 for £ > 1. 
This implies in particular that the DC (linear response) 
conductivity at non-zero temperature behaves as cr(T) ~ 
At 2 T 2A ~ 3 . These results obtain for the pure ring model, 
but will doubtless be altered somewhat in the presence 
of impurity scattering. This we leave for future study. 



V. INSTABILITIES OF THE BOSE METAL 

There are two classes of perturbations that one can add 
to the Bose metal fixed point (described by the quadratic 
action Sq) that can potentially destabilize the phase. The 
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first are terms involving the hopping of bosons. When 
relevant, such boson hopping terms "stiffen" the phase 
fluctuations of the bosons, leading to off-diagonal long 
range order and a superconducting state. As we shall see, 
the perturbative relevance of the boson hopping terms is 
determined by the Bose liquid parameters which enter 
into the fixed point action which characterizes the Bose 
metal phase. There are regions of parameters where all 
such boson hopping terms are irrelevant and the Bose 
metal phase is stable to such superconducting perturba- 
tions, as we detail in Subsection A below. 

The second class of perturbations involve hopping or 
motion of vortices, conveniently expressed in the dual 
representation. When relevant, these perturbations sig- 
nal a condensation of vortices which typically leads to a 
breaking of translational symmetry and drives the system 
into an incompressible insulating state. The presence of 
these instabilities and the precise form of the transla- 
tional symmetry breaking depend sensitively on the bo- 
son density, generally requiring boson densities commen- 
surate with the underlying lattice. In subsection B below 
we focus on half-filling (n = 1/2), and study the nature 
of the resulting commensurate insulating states. 

If the boson density is incommensurate with the lat- 
tice, on the other hand, small vortex hopping terms are 
unimportant. Provided the Bose liquid parameters are in 
the regime where Boson hopping is likewise irrelevant, the 
Bose metal exists as a completely stable critical phase. 
The gapless Bose metal is then a 2d analog of the stable 
Id Luttingcr liquid. 



A. Boson hopping and superconductivity 

Consider then the stability of the Bose metal phase 
in the presence of boson hopping operators. To be gen- 
eral, we consider processes where q bosons hop along a 
displacement vector s: 

£ t = -J252t s q cos[q(<p T -<p r+s )}. (118) 

r q,s 

To assess the perturbative relevance of such processes, 
we compute the two-point function of the tunneling op- 
erators, 

T 9 s (r) = cos[q(<p r -<p r+s )], (119) 

using the Gaussian action for the Bose metal. The 
row/column symmetries in the Bose metal phase greatly 
constrain the spatial dependence of these correlators. 
We consider first the special class of tunneling operators 
which hop bosons along the x or y axes with s = (to, 0) or 
s = (0, to). For this class of operators one finds a power 
law decay both in time and in one spatial dimension: 



with scaling dimension, 

f_ r dky l u(0,k y )sm 2 (mky/2) _ q*_ 
* qm ~2nJ_„ 27rV JC(0,k y ) \sm(k y /2)\ " 

Notice that these scaling dimensions vary as q 2 , 
and will generally increase (slowly A ?m ~ 
[q 2 / (27r 2 )] y/Uo //Co In |to|) with hopping distance to. As 
expected these scaling dimensions increase with U/K. 

For the remaining boson tunneling operators with vec- 
tor s connecting two different rows and two different 
columns, the correlator is further constrained by the 
row/column symmetry, being spatially local, 

<T g s (r,T)T g s (0,0))o = G(r)^o fors^O, (122) 

with G(r) ~ \t\- 2A «. 

General renormalization group reasoning implies that 
operators are irrelevant about a given fixed point when 
the associated scaling dimension exceeds the space-time 
dimension, A > D = d + 1, with dynamical exponent 
2 = 1. But due to the constrained form of the above 
correlators which only exhibit a power law decay in a 
reduced set of space-time dimensions, D rec i, one expects 
that the condition for irrelevance should be modified to 
be A > D re d. Thus, when A g)jn > 2 and A^ > 1, one 
expects that all of the boson hopping operators should be 
unimportant as one scales down in energy, and the Bose 
metal phase will be stable. This can always be achieved 
(in principle) by increasing the ratio of U/K . 

It is tempting to strengthen this expectation by con- 
structing an explicit renormalization group transforma- 
tion, but this is somewhat problematic due to the pe- 
culiar existence of zero energy states at both k x = 
and k y — in the Bose metal phase. One could try 
to integrate out gapped modes away from the zero en- 
ergy "cross" in the Brillouin zone, and then successively 
integrate out "shells" of modes pinching down onto the 
cross. A difficulty arises, however, in the rescaling trans- 
formation of the momenta, because the range of both k x 
and k y - the interval [— n, ir] - is not invariant. Similar 
but less severe difficulties are encountered when one tries 
to implement a momentum shell RG procedure for a 2d 
Fermi surface, due to the possible modifications of the 
shape/size of the Fermi surface. It might be possible to 
cirpjmvent this difficulty along the lines of Shankar et. 
alEEl, or by an RG procedure in frequency spaceH Some 
insight can be gleaned by ignoring the zero modes along 
the k x = axis, and considering a Id RG transforma- 
tion where the integration is over a shell of momenta in 
k x for all ky G [—it, 7r] and frequency, and then rescaling 
both k x and to (ie. with z = 1) but not k y . The result- 
ing perturbative (linearized) RG equation for the Boson 
hopping amplitude in the y— direction, t q ^ m = i™ y is then 
simply, 



<T™*(r, r)T q m * (O,0)) o = 5 yfi [x 2 + v 2 T 2 ]- A ^, (120) 



(123) 
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with z = 1 and A gjrra the scaling dimension given explic- 
itly above. This argument indeed supports the expecta- 
tion that weak Boson hopping will be irrelevant provided 
A 9 , m > D re d = 2. In the absence of a fully controlled 
2d RG procedure, we can verify that this conclusion in 
fact correct by resorting to perturbation theory In Ap- 
pendix ^], we do just this by formally computing correc- 
tions to the two-point function, (e lVr< - T ' > e~ liPr ( )) pertur- 
batively in powers of the Boson hopping t q —i <m —\. Sta- 
bility of the Bose metal phase requires that the long-time 
behavior be unmodified from that calculated within the 
Gaussian theory. Carrying out this expansion to second 
order in ti i, we show explicitly that this is in fact the 
case provided Ai,i > 2, thereby confirming the expected 
result. 

When the scaling dimension of a single boson hopping 
term is sufficiently small, Ai m < 2, one expects the Bose 
metal phase to be unstable to a superfluid state in which 
the bosons condense, with (e 1 ^) ^ 0.E3 If the original 
lattice Bosons are supposed to represent the Cooper pairs 
(say in a model of the cuprates) , then this will of course 
be the superconducting phase. 

B. Vortex hopping and insulating states 

1. Stability 

We next consider the effects of the various non-linear 
interactions involving the vortex operators in Eq. (|3§| ) 
and (|39|), which can potentially destabilize the Bose 
metal phase. At a generic density, for which n is irra- 
tional, all these cosine terms are oscillatory, and, if weak, 
cannot lead to any long- wavelength divergences. Hence 
we expect that away from commensurate densities, and 
provided A gTO > 2 for all q, to so that boson tunneling 
is irrelevant, the Bose metal is a stable zero temperature 
phase of matter. 

For rational densities, some of the vortex operators will 
be non-oscillatory, and must be considered more care- 
fully. In particular, for very commensurate boson den- 
sities (i.e. fi equal to a small-denominator rational frac- 
tion), it is not apriori obvious whether the Bose metal 
can even in principle be stable to both boson hopping and 
vortex operators. To demonstrate the issues, we present 
a stability analysis for the special set of rational fillings, 
n = 1/z, where z is a positive integer. 

At h = 1/z, the on-site terms take the form, 

= -^rf> 2g cos(2qtf+^:q/). (124) 

r q— 1 

Note that although for q =/= jz, with integer j > 1, 
these terms are spatially-varying, they are not wholly 
oscillatory, in the sense that they all have a non-zero 
spatial average for constant e.g. for prime z, and 
q =/= jz, [cos(^-xy)] xy — 1/z. Hence even for q =/= jz, 
they cannot be argued away simply on the grounds that 



they are oscillatory. Instead, to assess the importance 
of these perturbations, we again consider the two-point 
function of these operators evaluated with the Gaussian 
Bose metal action, 

Gf(r, t) = ( e *«*rM e -a«*>(0)) = s r , F(r), (125) 

With F(t) ~ e -9 2 VW^ln 2 (r) &% lajge t j mes _ Due to tue 

dual row/column symmetry this correlator is spatially 
local, and is also "short-ranged" in time, vanishing faster 
than any power law. The associated scaling dimension is 
thus infinite, and these operators are strongly irrelevant 
and will not destabilize the Bose metal. However, as we 
shall see, they will play an important role in determining 
which of the various insulating phases is selected when 
the Bose metal is driven unstable by the two-site terms. 

When the two-site terms in Eq. (^) are weak, we can 
drop those that are spatially oscillating and focus on the 
rest which take the form of vortex hopping and creation 
terms: 

C w = - w a ,±O a<± (r), (126) 

r,a,± 

with operators, a ,±(r) = coa(2q-d r ±2q'i3 r+s ), where "a" 
labels the various values of integers q, q' and the hopping 
vector s. Again, to establish the perturbative relevance of 
such terms in the Bose metal phase, we evaluate the two- 
point correlators with the Gaussian fixed point action. 
For the operators O+, the dual row/column symmetry 
is especially restrictive and we find that the associated 
correlator is once again spatially local: 

(O a ,+ (r, r)O o , + (0, 0)) = F(r)5 r , , (127) 

and "short-ranged" in time with F(t) ~ e -c a , q in (r) 
(with constant c a , q ), so that like the on-site terms, these 
operators cannot destabilize the Bose metal. Similar be- 
havior is found for C_, except for the special class of 
operators with q — q' , which correspond physically to 2q 
vortices hopping along a vector s. Of these, except for 
the special case with s = (to, 0) and s = (0, n), the two- 
point function is again spatially local, although it is now 
a power law in time, F(t) ~ |t| _2A . While such oper- 
ators could potentially destabilize the Bose metal if the 
power A < 1, they will generally be less singular than the 
remaining class of operators with vortices hopping along 
the x or y axes: 

C w = - W 2g , m 2J X! { c Os[2q{i!) X y ~ tfx+m.j/)] 
xy qm—jz 

+ cos[2g(0 BV -0 a , v+TO )]}. (128) 

The operators with qm ^ jz (with integer j) are truly 
oscillatory (i.e. have zero spatial average for constant 
i?) at fi = 1/z and have been dropped. We expect the 
coefficients W^qm to be positive (with the minus sign out 
front), since these operators will be generated at second 
order from the irrelevant on-site terms: w^q.m ~ v\ q . 
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For this last class of operators the two-point function 
is delta-correlated in one-spatial direction but a power 
law in the other and in time: 

(O y 2q . m (r,T)0%, m (0,0))o ~ S yfi (x 2 +«V)- A ^, 

(129) 

where 02 q>m (x,y) = cos[2q(tf xy - i? x ,y+m)]< The associ- 
ated scaling dimension is given in terms of the Bose liquid 
parameters: 



/C(0,fc,)sin 2 (mfc,/2) _ 
dky yU(0,k y ) sin(V2) - 2q ^ m) - 

(130) 

As for the boson hopping operators, the vortex hop- 
ping scaling dimensions vary as q 2 , and also increase 
slowly with increasing hopping distance, to. In general, 
we expect that for large z the most relevant of these 
will be it)2, zi with scaling dimension A 2 z = 2r) v (z) ~ 
2^/K.o/Uq In \ z\, so for sufficiently large z this can be ren- 
dered arbitrarily large, and hence all the vortex hopping 
terms can be made irrelevant for any W(k), /C(k). Thus 
for sufficiently large z, such that A 2qm > D rec i = 2 for all 
qm = jz, there is certainly a domain of stability for the 
Bose metal phase (i.e. where the Bose liquid parameters 
are further tuned to make all hopping terms irrelevant). 

For small z, however, this is not clear. Indeed, it 
is straightforward to show that for z = 2 (i.e. half- 
filling, n = 1/2), a choice of uniform functions U(k) — 
Wo,/C(k) = /Co does not lead to a stable regime. In par- 
ticular, in this case, the leading vortex hopping operator 
has scaling dimension A 2 2 = 2^(2) = A^JCq/Uq, while 
the leading boson hopping operator has scaling dimen- 
sion An = ^\/Uq/IQ). Hence A^Au = 4/7r 2 , which 
implies Min (An, A 22 ) < 2/ir < 2, so that at least one or 
the other operator is relevant. We have not determined 
whether a stable Bose metal phase might be possible at 
half filling when W(k),/C(k) are momentum-dependent. 



2. Instabilities 

As indicated above, for small z, either a vortex or bo- 
son hopping instability may be inevitable. In any case, 
it is interesting to study the nature of the state result- 
ing from relevant vortex hopping terms. Here, we briefly 
study the nature of the resulting phase for half-filling, 
taking into account the presence of the two most po- 
tentially relevant operators, with q = l,m = 2 and 
q = 2, to = 1. In general, at half-filling, there are Bose 
liquid parameters for which both w 22 and t± t i are rele- 
vant, and more complex behavior may well occur in this 
regime. We will however neglect boson hopping com- 
pletely, as appropriate for large U/K,. 

To this end, let us assume that u> 2;2 is the most rele- 
vant operator, with A 2 2 < 2. We also assume that 
is the next most relevant (it could be irrelevant, but still 
the least irrelevant remaining operator). Then provided 
all bare couplings are small, we imagine integrating down 



in energy until w 2t2 becomes comparable to the Gaussian 
terms in So- This requires k; 2;2 A A2 . 2 ~ W(7r,0)A 2 , which 
always occurs for sufficiently small A with A 2 2 < 2. At 
this point it seems appropriate to minimize the potential 
-w 2 ,2[cos(2'd X y - 2-d x+2 ,y) + cos(2'd X y - 2-d x , y+2 % simul- 
taneously with Sq. This minimization somewhat under- 
estimates fluctuation effects, which will be commented 
upon later. The most general form for $ xy which mini- 
mizes the W2.2 term and keeps the Gaussian action small 
can be written 

v r = e (r) + (-rre^r) + (-l)»e 3 (r) + ^ + ^, 

. (131) 

where ©o/i/2( r ) are slowly-varying functions of space. 
Inserting this expression for d into the action including 
So and the u> 2)2 term only gives 



<->o 



•err 



7T 2 /C 



(drOo) 2 



TT 2 ^ 



[(9 r 6i) 2 + (d T e 2 ) 2 



^(ve o ) 2 + ^(9 x 0i) 2 
^(d x e 2 ) 2 + ^(d y Q 2 ) 2 



(132) 



where JC V = /C(tt,0), A+ = 8U{n, 0)/tt 2 , A. 



_ u> 22 A 



w. 



2/(2-A£, 2 ) 
2,2 



-A" /(2-A" ) 



2 - 2 > where A 



is the reduced low energy momentum cutoff. Note that 
we are not rescaling any fields or coordinates in this 
schematic RG treatment, just taking into account the 
fluctuation corrections to u> 2j2 . Thus the low-energy con- 
tinuous fluctuations around the minima are described by 
three massless fields Qj (j = 0,1,2), which unlike $ in 
the absence of w 2i2 have an ordinary "relativistic" dis- 
persion. In addition to these continuous variables, the 
discrete degeneracy of distinct minima allowed by peri- 
odicity of the cosine is indexed by the integers a, b. 
Having taken into account w 2 , 2 already in Eq. (132), 



we next include the effects of a weak, renormalized 104,1 
coupling, which after renormalization is of order ~ 

i«4,iA a m - W4,i(w 2 , 2 M7r,0)) A 4.i /(2 ~ A2; . 2) . If all the 
bare non-linear couplings are small, then after renormal- 
ization, this will be the largest remaining term, which 
is why we treat it next. Inserting the decomposition of 
Eq. (|13ip into the four vortex hopping term, one obtains 



•off 



-<1 



[cos89i + cos80 2 ] 



(133) 



assu ming small gradient s of @i 



as mandated by 

Eq. (11321). Eqs. (Il32l,|l33j) describe a 3d "height" model 



/2= 



for ©i, 8 2 . As is well known in such models, the fluc- 
tuations of the free scalar fi elds ©1, 9 2 are bounded, so 
that the cosine terms in Eq. 133 are always relevant and 
pin these fields to integer multiples of 7r/4, which we can 
take, without loss of generality, to be zero. Thus the net 
effect of W4.1 is to leave only ©0 as a low-energy mode. 
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Setting therefore 0i 
phase, we have 



9 2 = in the pinned ("flat") 



(134) 



Here, since the fluctuations of 0o are likewise bounded, 
we have replaced it by its average (zero mode) value 
@o — * i^o- This average can be de term ined by mini- 
mizing the on-site Lag rang ian in Eq. ( 124 ). In principle, 
the parameters in Eq. (124) should also be renormalized, 
e-g- v 2q — > v 2q ~ U2 9 exp[-g 2 c g ln 2 (w 2 ,2/^(7r,0))]. This 
gives, 



^2ii9 



-1) 



ah 



(135) 



for V4 > —V2/8 (with vi > assumed) and 



° = (-l) ab [(v 2 /8\v A \) ± iVH^F], (136) 

for U4 < —1)2/8. 

The spatial orderings that are implied by these mean 
field solutions follow readily from the expressions in 
Eq. (^,|l^) relating the bare boson density and plaque- 
tte energy density to the "low-energy" field '&. In the 
former case above, the boson density is uniform since 
sin(2q'& Xi y) = 0, whereas there is a plaquette energy den- 
sity wave with, 



K-i) 



(x+a)(y+b) 



(137) 



The four states with a, b = 0, 1 correspond to the four 
plaquette density wave states of the original lattice Bo- 
son model in which one out of every four plaquettes is 
resonating more strongly. 

When Vi < —1)2/8, on the other hand, since sin(2$) ^ 
the plaquette energy density wave states (with reduced 
order c\ — > (£>2/8|w4|)c|) are co-existing with a charge- 
density- wave (CDW) state ordered at (7r,7r): 



5n x+hy+ , ~ ±cWl-(v 2 /8v4) 2 (-l) x+y . (138) 

In the limit in which £4 — > —00, the amplitude of the pla- 
quette density wave vanishes, leaving only CDW order. 
We believe that an absence of a pure CDW state more 
generally is an artifact of the mean-field treatment being 
employed, which ignores all fluctuations in the i9 field. 
In particular, one can imagine domain walls forming be- 
tween domains of different plaquette-density-wave order, 
which cost an energy E wa u ~ u^^fe/'tU) 2 for small w 2 ,2 
and large \vi\. When this energy is small, fluctuations 
will undoubtedly disorder the plaquette-density-wave or- 
der, leaving the pure CDW. 

Our quantum Monte Carlo simulations on the boson 
ring model in the dual representation reveal that in the 
large U limit with U/K > (U/K) c w 2.5, the Bose metal 
phase is unstable to the formation of a CDW state. This 
is apparent in Fig. ^ which shows that the density struc- 
ture factor peak at (tt,tt) grows as L 2 , indicating the 
presence of long-ranged CDW order. This order is, how- 
ever, weak, with a staggered magnetization to ~ 0.1 (see 
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FIG. 6: Finite size scaling of the density structure factor peak 
at (it, 7r), well within the incompressible phase. The staggered 
magnetization m = (m 2 ) 1 ^ 2 clearly saturates to a finite value 
1/L — > 0, indicating a non-zero (tt, tt) charge density wave 
order parameter. 
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FIG. 7: Vortex correlation function G$(0,t) in the incom- 
pressible CDW phase. The correlation function saturates as 
seen from the finite size scaling plot in the inset indicating 
that the vortices condense, leading to a conventional insulat- 
ing state. 



Fig. |^)), much smaller than the classical value to = 0.5. It 
is also instructive to examine the vortex two-point func- 
tion in the CDW state, which is shown in Fig. [?]. Notice 
that the vortex correlation function is long-ranged indi- 
cating a "vortex condensation" ((e ,e ) 7^ 0), as expected 
in such a conventional insulating state. 

Both the plaquette and charge density wave states will 
be insulators, with a charge gap. This follows since in 
both cases the field ■& is "pinned" by the cosine poten- 
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tials and is not fluctuating. Since the bare boson density 
is ri = A xy 9/n, adding a particle at the origin can be 
achieved by shifting 9 xy — > 9 xy + tt for all x, y > 0. This 
will cost a finite amount of energy (coming from the pla- 
quette at the origin) when the field $ is pinned. 



VI. EXACT WAVEFUNCTION 

In this Section we consider a variant of the boson ring 
Hamiltonian which allows us to obtain an exact zero 
energy wavefunction when the couplings are carefully 
tuned. We then perturb away from the soluble point, us- 
ing the exact wavefunction to compute properties of the 
adjacent quantum phase. Specifically, we find a transla- 
tionally invariant fluid phase with a finite compressibility, 
behavior consistent with both a superfluid and the Bose 
metal. But a calculation of the boson tunneling gap in 
a finite system shows a 1/L scaling as expected for the 
Bose metal, and inconsistent with the 1/L 2 dependence 
of a 2d superfluid. Thus, we confidently conclude that 
the Bose metal phase exists over a finite portion of the 
phase diagram adjacent to the soluble point. 

The model we consider is a "hard-core" version of the 
boson ring model, in which only zero or one boson is 
allowed per site, tfb = 0,1. In the usual way we represent 
the hard-core bosons as Pauli matrices, 



b 

n = tfb 



~ 2 {l + °% 



(139) 
(140) 

(141) 



where a is the standard vector of Pauli matrices, and 
a 



- \{p x ±ia v ). The Hamiltonian we consider consists 



of two terms: 



H i/2 - Ej - ■h(cr+ y a x+l y a^ +l y+1 a x y+1 + h.< 

xy 



+u i Pm P (x,y) 



(142) 



where 



1 

Tg 



e n 11 

a=±l /3,7=±1 



(143) 

The first term with J4 > in Eq. (142) is the hard-core 
analog of the ring hopping term proportional to K in the 
rotor model, Eq. ([!]). The operator Pfn p (x,y) is a pro- 
jection operator onto the two flippable configurations of 
the square plaquette whose lower-left corner is at the site 
(x, y). For 7/4 > 0, this term competes with the ring term 
by disfavoring configurations with flippable plaquettes. 

Remarkably, the ground state of ifi/2 can be found ex- 
actly for the special Rokhsar-Kivelson (RK) point U4 — 
J4, following a general construction in the spirit of the 
Rokhsar-Kivelson point of the square lattice quantum 



dimer model^l and employed more recently in Ref. [l6| 
for a similar spin model on the Kagome lattice. The so- 
lution can be seen by rewriting if 1/2 as 

#1/2 = E^^H* 7 ^ 1 - II °Z+v,v+v) - v }' 

xy u,v— 0,1 

. - r (144) 

where v = — measures the deviation from RK point. 
For v — 0, an obvious (zero energy) ground state of if 1/2 
is the fully-polarized state 



io)=ni< 



=1). 



(145) 



A useful alternate representation follows by rewriting 
\<j x = 1) = -^(W z = 1) + \a z = —1)), and expanding 
out the direct product, 



|0) 



1 



2 N/2 



E lite 

{cr XH =±l} xy 



(146) 



which demonstrates that |0) is a uniform real superposi- 
tion of all configurations in the S z (boson number) basis. 

This state, however, has uncertain a z (boson number), 
and so can be decomposed into many distinct ground 
states by projection. In particular, at a given average 



■ryl 



2(n x 



1 = m, we may project onto spatially 



uniform states according to 

L 



|0;m) = ^= f[ P x (m)Py{m)\0), (147) 

where Z is a normalization constant, since the 
row/column projection operators, 



P x (m) 
P y (m) 



2 IT 



d6 



, 2tt 
2tt 



exp 



exp 



E «u - m ) 

v 

-i<t> y {tey - " 



,(148) 



commute with ifi/2 and amongst themselves, reflecting 
the conservation of boson number on each row and col- 
umn by the ring dynamics. 

Of course, other non-uniform ground states may be 
obtained by choosing m differently on different rows and 
columns. This vast degeneracy signals a pathology of 
the RK point, which is in fact at the boundary of a first- 
order transition to a phase-separated "frozen" regime, 
for v < 0. We therefore focus instead on the behavior for 
infinitesimal v > 0, which splits this degeneracy. For bo- 
son densities near half-filling, i.e. \m\ <ti 1/2, we expect 
that the uniform states will be favored energetically as v 
is increased to slightly positive values, since these states 
have in this density regime more flippable plaquettes (see 
below). For larger \m\, this assumption is certainly vio- 
lated, however, since e.g. at very low boson density the 
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ring moves clearly do not connect all possible configura- 
tions. At m = 0, however, the set of uniform configura- 
tions does form a single ergodic component under the ring 
move (as can be straightforwardly shown numerically^ 
and probably argued analytically). So we expect the set 
of uniform states to be an adequate description for small 
\m\. 

To make contact with the Bose metal fixed point de- 
scription, we calculate some energetic properties for in- 
finitesimal positive v using first-order perturbation the- 
ory. From the splitting of the different projected states at 
0(v), we calculate two quantities: (1) the ground state 
energy density e(m) = E(m)/L 2 as a function of bo- 
son density n = (m + l)/2, from which we obtain the 
compressibility n(m), and (2) the "single-particle" gap 
Ai(m, L) for a finite-size L x L system, essentially the 
addition energy for a boson onto a particular row and 
column in the grand canonical ensemble, defined more 
precisely below. 

Consider first the ground state energy density. For the 
projected state at magnetization m, the first order energy 
shift is 



E(m) = -u^(0;m|P fli p(x,y)|0; 



xy 

L 

= -|^(0|P fl i P (z,y) J] 4(m)4(m)|0),(149) 

xy x,y—\ 

where the latter equality follows from the fact that Pfli p 
is diagonal in the a z basis and he nce c ommutes with the 
projection operators. Using Eq. (145), one readily sees 



that the energy shift can be rewritten in a form reminis- 
cent of classical statistical mechanics, 



E(m) 



'^(Pfiip(^,y;{o-})) {c 



(150) 



where 



E 

{<j J . y =±l} x,y=l 



O H P x (m;{a})P v (m;{a}). 



(151) 

Here -Pflip, P x and P y are the classical functions obtained 
from Pflip, P x , P y , respectively, by replacing the operator 



Vxy by 

and Z is chosen such that (l){ CTx 



,} 



1. 



We see that Eq. ( 151 ) simply defines an expectation 
value for an "infinite temperature" Ising model with a 
constrained magnetization on each row and column. By 
proceeding along similar lines, one can calculate the en- 
ergy shift for a state with an extra boson on a single 
particular row and column, and hence obtain Ai. Amus- 
ingly, both problems can be solved exactly using a saddle- 
point technique (see Appendix ^) . The results demon- 
strate that the constraints are "nearly irrelevant" , and 
the lattice gas behaves nearly as its unconstrained coun- 
terpart in the large system limit. In particular, the en- 
ergy density, as L — ► oo, becomes 



(m) 



2\2 



(152) 



This can be easily understood by assuming that each site 
is completely independent, and that the only effect of the 
constraint is to determine the relative probabilities of the 
two spin states, to wit = (1 + m)/2,p| = (1 — m)/2. 
On a given four-site plaquette, only the two (of 16 total) 
configurations in which the spin alternates around the 
plaquette are flippable. Hence the average flippability 
per plaquette is 2p|p| = (1 — m 2 ) 2 /8, in agreement with 
the result above. 

The "single-particle" gap is somewhat less intuitive. 
To be precise and avoid ambiguities due to an unspecified 
chemical potential, we define 

1 



A x (m) = -[E 1 (m; 



-l)-2P 1 (m;0) + Pi(m,-l)], 

(153) 

where E\ (m; A) is the ground-state energy of a state in 
which the number of bosons on row and column x = y = 
1 is increased by A relative to the the uniform state at 
magnetization m: 

|0;to,A) = -=P x = 1 (m + 2\/L)P y=1 (m + 2\/L) 
v 



x II Px(m)P y (m)\0), 

x,y=2 



(154) 



and Pi(m; A) 
finds 



~ v Y^x V (ty m ' A|Pfl ip (a;, y)|0; m, A). One 



Ax(m) = ^(l 



m 2 ) + 0(l/L 2 ). 



(155) 



These two quantities can be compared to the gen- 
eral predictions expected from our theory of the Bose 
meta l ph ase. Consider first the compressibility. From 
Eq. (152), one has fi = —de/dn = —2de/dm — vm(l — 
m 2 ), and hence 



dn 



1 



1 



2v 1 — 3m 2 



(156) 



Note that the compressibility diverges and becomes neg- 
ative for l/\/3 < \m\ < 1, indicative of an instability 
of the uniform state well away from half-filling. For 
|m| < l/VS, however, k is finite, consistent with the 
Bose metal phase. Indeed, from our general Harmonic 
description of the Bose metal phase as in Eq. ??, we 
have k- 1 = U{k = 0). 

In order to rule out a superfluid phase, we now com- 
pute the single-particle gap, Ai(L), using the Bose metal 
Gaussian fixed point theory to show that it varies as l/L 
as in Eq. (155). This gap can be extracted from the spa- 
tially local correlator, 



— /pWr(r) -ivr(0)\ 



(157) 



evaluated in a finite L x L system. From the spectral 
representation, one has, at zero temperature, 



Ai(L) = lim - 



In G v (t) 



(158) 
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Performing the Gaussian integral using the effective ac- 
tion in Eq. (|43|), one has 



In G v (t) 



du 1 \ - 



(1 



(159) 

where the wavevector sum is over inequivalent values 
with k x ,k y G (2ir/L)Z in the first Brillouin zone, and 
with the mode energy E^ given explicitly in Eq. ( p7| ) . As 
t — > oo j the oj integral in Eq. (159) is dominated by those 
terms for which E^ = 0. This occurs along the Bose sur- 
face, i.e. for k x = or k y = 0. In the finite size system 
there are 2L — 1 w 2L such points (for L> 1), hence 



2U 



where 
U = 



71 dk 
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-U{k x ,0) 



dk 
.Ytt 



^W(0, ky 



(160) 



(161) 



is the Bose-surface average of U(k). Thus, we finally 
have, 



lnG (O,r) 



L ' 



(162) 



giving the l/L behavior for the single-particle gap: 



Ai 



U 
L' 



(163) 



In a 2d superfluid, the single-particle gap is much smaller, 
vanishing as l/L 2 in the large L limit. 

The agreement of the finite compressibility and l/L 
scaling of Ai between the soluble model and our fixed 
point theory of the Bose metal strongly argues that the 
ground state of H 1 / 2 is the Bose liquid for < v <C 1 
and |m| < l/VS. If this postulate is correct, certain 
combinations of the Bose liquid function U(/k) can be 
obtained explicitly for the soluble model. In particular, 
we find 



W(k = 0) = 2v(l-3m 2 ), 
U = 2v{l-m 2 ). 



VII. DISCUSSION 



(164) 
(165) 



A. Fractionalization, the Z2 gauge theory and the 
high-r c cuprates 

In this paper we have described a novel phase of quan- 
tum matter, the Bose metal, which we argue occurs in a 
class of square lattice boson ring models with XY symme- 
try. Having done so, it is reasonable to reflect upon the 
context in which such models are physically appropriate 
and the consequences perhaps observable. As discussed 



briefly in the introduction, the primary motivation for 
these models comes from the high temperature cuprate 
superconductors. Both the remarkably high critical tem- 
peratures and the strange behaviors in the "normal" state 
of these materials motivated a number of theorists early 
on to the radical suggestion that soifccbacge separation 
might underly these peculiarities P^P^P^P 8 ! In this pic- 
ture, the electron charge, liberated both from its spin and 
its Fermi statistics as a "chargon" (or "antiholon"), is rel- 
atively free to Bose condense to form a superconducting 
state. At higher temperatures, or when the material is 
under-doped, the chargons and spin-carrying fermionic 
"spinons" form a novel fluid, which would no doubt be- 
have very differently from a normal metal. 

Subsequent theoretical and consequent experimental 
work has both advanced the status of these ideas and 
posed a severe challenge to theip-, applicability to the 
cuprates. A recent formulationtML3 of interacting elec- 
trons in terms of spinons and chargons minimally cou- 
pled to a Z2 gauge field has helped to support the con- 
cept of electron fractionalization with a very concrete for- 
mal framework. The essential elements of the Z2 gauge 
theory are reviewed in Appendix As with., earlier 
U(l) and SU{2) gauge theory formulationsOE3cil, the 
"spinons" are taken as Fermions carrying the spin of 
the electron but are electrically neutral. The bosonic 
"chargons" carry the electrons charge. The Z2 gauge 
theory provides a convenient phenomenology for describ- 
ing a fractionalized phase in which the spinons and char- 
gons are deconfined, and live as well-defined particle ex- 
citations. Because of its concrete formulation, and the 
virtue that the Z2 gauge theory has a well-understood 
confincmcnt-dcconfincmcnt phase transition, it is possi- 
ble to calculate from it in a simple way the qualitative 
(universal) properties of the fractionalized state. One 
of the most fundamental of these properties is the ex- 
istence of gapped topological excitationsQaiJ called "vi- 
sons" . The visons act as sources of Z2 flux, and may 
be thought of as remnants of unpaired superconducting 



which in turn is 
Understanding 



vortices in the fractionalized insulator', 
viewed as a paired- vortex condensates . 
these excitations led to a rather direct proposal for a "vi- 
son trapping" experiment ,E9 designed to trap and detect 
visons by cycling through the supjpconducting-normal 
transition. Current experimentst£IE3 apparently imply 
that the gap to such vison excitations in the insulator, if 
it is non-zero at all, is less than 190-ftT. Since the presence 
of a vison gap is a necessary condition for the existence of 
a true spin-charge separated ground state, this unnatu- 
rally low energy scale presents a difficult obstacle to the- 
ories of fractionalization as applied to at least these par- 
ticular cuprate materials (samples of underdoped BSSCO 
and YBCO). 

The spin-charge separation scenario is nonetheless ex- 
tremely appealing theoretically, and it is interesting to 
consider the possibility of retaining some degree of this 
physics locally. One may imagine for simplicity an un- 
doped model which interpolates between at one extreme 
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a conventional "Hubbard-like" insulator with gapped 
charge degrees of freedom and Heisenberg-interacting 
spins, and at the other extreme a fully fractionalized in- 
sulator. The Z2 gauge theory is well-suited for this pur- 
pose. The interpolation between the above two limits is 
accomplished in this model by varying a coefficient K, 
which controls the strength of fluctuations in the gauge 
field. The deconfined phase is obtained in the large K 
limit, while we consider here the first deviations away 
from K = 0, which is deep within the confined, Hubbard- 
like phase. As we show in Appendix when the gauge 
theory is deep within it's confined phase (with K small), 
the gauge fields can be formally integrated out and one 
recovers a Hamiltonian expressible in terms of electron 
operators and composites built from the electron such as 
the spin operator, together with a Cooper pair field. In 
the spin sector, we find that the leading terms obtained 
in this limit are an antifcrromagnetic nearest-neighbor 
Heisenberg exchange (for K = 0) and a plaquette ring 
term, 

i^=J5^ s D , (166) 

□ 

with J* - O(K) (J» = K(t 4 s +A 4 )/(h + U) 4 , in terms of 
the parameters in Appendix Here, lZ s a is a plaquette 
ring operator defined in terms of the four spins on a given 
plaquetteBEi 

4 

W n = S i -S i +4(Si-S a )(S3-S 4 ) 

«3'=1 

+ 4(Si • S 4 )(S 2 ■ S 3 ) - 4(Si • S 3 )(S 2 • S 4 )(167) 
In the charge sector, to O(K) one finds the leading term 

H^ = -J^TZ c n , (168) 

□ 

with - 0{K) (Jg = Ktj/(h + [/) 4 , derived in Ap- 
pendix |F^). Here, TZ^ is a ring operator in the charge 
sector expressed simply in terms of Cooper pair opera- 
tors, e l * r = b"l as, 

4 

~R-n = cos(^i -<j)j) +cos(A a ,j / r ). (169) 

i<j=l 

Notice that the second term above is precisely of the form 
taken in our starting Hamiltonian. As K is increased 
from zero, such ring-exchange processes play more and 
more important roles. We thus view the ring Hamiltonian 
studied in this paper as a suitable model that offers an 
intermediate ground between the fully spin-charge sepa- 
rated scenario and conventional phases of matter. 

Should this bosonic ring model for the charge sector 
have some relevance for the cuprates, what might be the 
consequences and interpretation? As we have seen in 
Sec. M, the ring model sustains the Bose metal phase for 



sufficiently repulsive interactions (large U) and incom- 
mensurate densities. Indeed, as we saw in Sec. [v], a sim- 
plistic estimate indicates that there is no stable regime 
for the Bose metal at half-filling. Hence, if we assume 
the strong interaction condition obtains, one would ex- 
pect an evolution from a conventional insulator with a 
charge gap at half-filling to the Bose metal upon dop- 
ing sufficiently away. If one further presumes (as seems 
natural) that the generalized "stiffness" K, increases with 
doping, then one would expect further doping to lead to 
a superconducting (boson hopping) instability. This sce- 
nario thus naturally associates the Bose metal with the 
pseudo-gap regime of the high-T c cuprates. 



B. Extensions 

To determine if there is any truth to the above scenario 
requires considerable extensions of the present work. 
Most significantly, spin and quasiparticle degrees of free- 
dom are important components of the high-T c materials, 
and should be incorporated into the description. It will 
be interesting to consider interactions of the Fermionic 
quasiparticles with the strongly fluctuating collective 
modes of the Bose metal. It seems to us quite likely that 
the Bose metal can remain stable in the presence of these 
interactions, while perhaps at the same time producing 
rather strong modifications of the fermionic degrees of 
freedom. In any case, it should be possible to consider 
photoemission spectra and local electron tunneling den- 
sity of states in a model with the quasiparticles coupled 
to the Bose metallic modes. 

We have also left a number of issues within the purely 
bosonic description unanswered. To understand trans- 
port measurements generally will require an understand- 
ing of the effects of disorder. Several experiments at- 
tempting to access the "normal" state of the cuprates 
employ large locaLar uniform magnetic fields to suppress 
superconductivity^. Such magnetic fields indeed tend 
to suppress the superconducting instabilities of the Bose 
metal, and one should explore the ramifications of a Bose 
metallic description for e.g. transport measurements in 
high fields. 

There are a number of other possible applications of 
the present formalism to frustrated magnets with appre- 
ciable ring exchanged. Magnetic ring exchange processes 
are believed important for instance in the^Wic'ner crystal 
phase of the two-dimensional electron gasBcJ at r s > 40. 
The effect of such processes are generally difficult to an- 
alyze, but a perhaps illuminating approach may be to 
consider the easy-plane limit of such spin ring models. A 
straightforward calculation shows that such a limit recov- 
ers XY ring models of the sort studied here (representing 
the spins by hard-core bosons). The plaquette duality 
constructed here for the square lattice can be straight- 
forwardly generalized to other lattices, such as the trian- 
gular and Kagome cases. For the bosonic Kagome case, 
one can directly in this way construct a dual formulation 
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from which it is st raig htforward to construct the visons 
and spinons of Rcf. Id Whether the spinons are confined 
or deconfined in the pure XY ring model on this lattice 
can be easily established numerically using the dual 9 
variables, and this investigation is undcrway.E3 On the 
triangular lattice with four-site (parallelogram) ring ex- 
change, the plaquette duality shows that there is no sta- 
ble metallic state in this case, but a numerical approach 
will be necessary to determine the ultimate nature of the 
ground state. 

On the purely theoretical side, we have introduced the 
mathematics necessary for a 2 + 1 dimensional "bosoniza- 
tion" scheme, applied here to "bosonize bosons" in terms 
of collective (p, & modes. We expect it should be possible 
to bosonize fermionic ring models in a similar way. Con- 
sider spinless fermions {/i,/J} = Sij (not at this point 
to be associated with cuprate quasiparticles) living on 
the sites of the original square lattice. One may formally 
define Jordan- Wigner hard-core bosons by introducing a 
"string" S winding around the lattice, 



fr — b v S r 



(170) 



where 



IK 



y>>y x' 



With this definition, the boson operators obey proper 
canonical commutation relations at different sites, and 
have a hard-core interaction on the same site. For an or- 
dinary dynamics, this extremely non-local string presents 
unsolvable difficulties for analytic treatment. However, 
using the exact plaquette duality appropriate for ring 
dynamics, one finds that one can write 



v+i 



(172) 



Here we have ignored any possible boundary terms. Up 
to this proviso, the string becomes local in the 9 vari- 
ables! Thus, in a rotor representation with b — e 1 ^ (this 
requires a large Un(n— 1) term to maintain the hard-core 
constraint), one has 



Clearly there are a range of applications to be explored 
based on the present work. We hope that some of these 
may ultimately enhance our understanding of experimen- 
tally accessible strongly correlated materials. 
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APPENDIX A: TROTTER DECOMPOSITION 

The usual Trotter decomposition of Z a gives 

T=0 



|0rr)) 



z = e n^+ei e£ * E ' cos(OTWr)e ~^ E ' [A * Ber " 

{e„}r=o 

. ( A1 ) 

where we have used e <C 1 to separate the imaginary- 
time evolution operator into to factors. The second can 
be directly evaluated to give 

Z= ]T JJe-^E^^- R l 2 n<<W|e eifcos ™|0 r 

{8 rT }T=0 r 

(A2) 

The latter matrix element can be written 



9 /| e eKco S (TJV r )|^ = (6>'|l + eXcos(7riV r )|6>) 

eK 



exp 



(A3) 



correct to O(eK). This leads directly to Eq. (El 



,v+i 
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which is directly analogous to the bosonization formula 
for spinless fermions in one dimension! This approach 
may lead to an understanding of novel non-Fermi liquid 
states in fermionic ring models. 



APPENDIX B: ASYMPTOTICS 

In this appendix we calculate the asymptotics of the 
correlator c xv (t, t±, T2), introduced in Sec. IV. Using the 
general rules of Gaussian theories, one has 



Cx, y (r,ri,T 2 ) = 2^(^(r)+i9oo(0)-^o(Ti)-^(r 2 )) 2 ) o 
= 27r 2 Re / £(k) i '- i — 



e -E k \T-r 2 \ \ e -ik x x_ ^ e --Ek|r 2 | 



e -B k |r-Ti|\ e -ik y y 

'—m 
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The above integral was considered earlier for the special case r = t\ — t 2 — 0, for which it was shown to grow 
logarithmically with x at fixed y (and vice versa). This growth with x and y persists when fiol r l; ^ol T i|i ^o| T 2| \%\ 
(or \y\). To see this, we differentiate 



d x c xy = 27r 2 Re 
JC{0,k y ) 



ik x e- tk * x )C(k) 
El 



g-S k |ri| _|_ e -E k \T-r 2 \ _ e -ife 



I 4_ e -Bk|ri-r 2 | 



w{k y ) 



x 2 +cj 2 (k y )rf + x 2 +u 2 (k y )(r-T 2 ) 2 C ° SkyV 



X 2 +LU 2 (k y )T 2 X 2 +LU 2 (k y ){T 1 -T 2 ) 2 



02) 



where uj(k y ) = Sl(0,k y )\2sm(k y /2)\, and in the second 
line we assumed (as valid for large |ar|) the k x integral is 
dominated by small k x and integrated it in that regime. 
Hence for \x\ » S7o|t|i 3 ^o| r l, one has 

c xy ~ <&] v (y)ki\x\ + f(T,n,T 2 , y), (B3) 

where / is an unknown (but ^-independent) function of y 
and the imaginary times. We expect that c xy should cross 
over to a In 2 |OoT"i | behavior if any of the fio|rj| become 
large compared to the spatial distances involved. How- 
ever, it is unclear without further calculation whether 
this crossover occurs when the imaginary time argument 
becomes large compared to the smaller (e.g. y above) or 
the larger (e.g. x above) of the two spatial coordinates. 
To ascertain this information, we differentiate with re- 
spect to the largest time (which we choose positive with- 
out loss of generality), e.g. for r > |ri|, \t 2 \, 



d T c xv = 2tt 2 f /C(k) 



-E k (T-Ti)-ik y y 



_|_ e --E k (T-T 2 )-ifc x z _ e -E k T-i(kxX+k y y) 



27T' 



(B4) 



where in last line we have kept the dominant term for 
\x\ ^> \y\,il T. If we assume \y\ ~ 0(1), then for Q (t — 
t\) S> 1, the integral is dominated by || w (due to 
the logarithmic singularity in the collective density of 
states), and one may approximate e~ %kyV w 1, leading to 
d T c xy ~ 2 v //C /^oln|»o(r - n)|/(r - n), whence / ~ 
y/ /Co ]Uq In 2 |Oo(r — n)|. More generally, one can show 
that this behavior obtains provided fio( T — r i) 3* \y\- By 
identical arguments, / also grows like v^Mln 2 \fl T 2 \ 

for Sl 1 T"2 1 3> \y\- Thus, the integral defining C X y (t) is 
dominated by t\ m t and r 2 w 0. Hence, we can focus 
on this case, and consider the behavior in general as a 
function of |x| S> 1 and t» 1, but with no particular 
relation between them. One finds 

/C(k) 



,(r, r,0) « 47r 2 Re 



(1- 



D -ik y y\ 



M_ e -E k \r\-ik x x\ 



/k #k 

It is straightforward to see that this integral recov- 
ers the result, Eq. ( [B3| ) for |x| ^S> |r|, and conversely, 
c xy (r,T,0) ~ 4^(y)lnr for |r| S> Putting all the 
above analysis together, one arrives at the general ap- 
proximate result of Eq. (p2|) in Sec. IV. 



APPENDIX C: RELEVANCE OF HOPPING 

In this Appendix we compute the corrections to the 
Bose metal imaginary-time equal-space correlator, 

G(r) = (e^'We^t )), (CI) 

to leading (quadratic) order in the Boson hopping term: 

S t = -ty / COS(lfi r + y ~ (f r ). (C2) 

Perturbing around the Gaussian Bose metal theory gives, 
G(r) = G (r) + G 2 (r) + 0(t 4 x ), with 



G 2 (r) = i( e ^M-^(°)5 2 ) c 



(C3) 



where the connected correlator is taken with respect to 
the Gaussian action. Using Wick's theorem this can be 
re-expressed as, G 2 (t) = f 2 Go(r)/(r), with 

I(r)= J2 I V {x 12 ,y 12 )[e- A -l], (C4) 
where the function A = A(xi , x 2 , y\ , ri, t 2 , t) is given by, 



A= e lkri - MTl B(k,w,r)(l-e 

J k,u; 



-ifexu+imj) + h.C, 

(C5) 



r iur _ - e iky ) 
S(k,a;,r) = W(k)^ 2 | l „ 2 - J . (C6) 

CJ -t- £, k 

Here, £12 = cci - x 2 , T i2 = n - t 2 and 2? (a:,T) ~ 
(x 2 +t 2 )- a is the two-point correlator for the tunnel- 
ing operator evaluated in the Bose metal. To isolate the 
t dependence of I(t), we re-scale x\,x 2 ,ti, t 2 by r, and 
k x and w by 1/t (leaving y and fcj, alone). For large r 
this gives, I(t) ~ t 2 ( 2 ~ a '/(1), or equivalently, 



^~t 2 /(l)T 2 ( 2 -*>. 

G (r) 



(C7) 



We thus conclude that provided 1(1) is finite, the large r 
behavior of the Bose metal correlator, Gq(t), is unmod- 
ified if the scaling dimension A > 2. To show that 1(1) 
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is finite, it suffices to expand the exponential to second 
order in A, which gives: 

id) = 8 / ^pfe.^ fi-^y-r^ , 

(C8) 

with 

V(k x ,u)= [ P (x,r)(l-e^™)). (C9) 

Provided A > 1 the integrals over x and r converge, 
givin g V (k Xl uj) ~ (fc 2 , + w 2 )^ 1 . Inserting this into 
Eqn. |C8j, one readily sees that the k and uj integrals are 
likewise convergent, confirming that 1(1) is finite when- 
ever A > 1. 



APPENDIX D: Z2 GAUGE THEORY AND RING 
EXCHANGE 

A recent formulation! of interacting electrons in two 
dimensions has been developed which re-expresses the 
electron operator in terms of "spinons" and "chargons" 
which are minimally coupled to a Z-i gauge field. As 
with earlier U(l) and SU(2) gauge theory formulations, 
the "spinons" are taken as Fermions carrying the spin 
of the electron but are electrically neutral. The bosonic 
"chargons" carry the electrons charge. The Z2 gauge the- 
ory provides a convenient phenomenology for describing 
a fractionalized phase in which the spinons and chargons 
are deconfined, and live as well-defined particle excita- 
tions. The deconfined phase is most readily accessed by 
increasing the strength of a term in the gauge theory 
Hamiltonian (with coefficient K, below) which suppresses 
the fluctuations in the gauge field. When the gauge the- 
ory is deep within it's confined phase (with K = 0), on 
the other hand, the gauge fields can be formally inte- 
grated out and one recovers a Hamiltonian expressible 
in terms of electron operators and composites built from 
the electron such as the spin operator, together with a 
Cooper pair field. Here, we show that upon integrat- 
ing out the gauge field with small but non-zero K, one 
generates ring-exchange terms in the effective Hamilto- 
nian. In the spin-sector of the theory, the leading order 
ring term involves four-spin operators around an elemen- 
tary (square) plaquette, and is the familiar quartic form 
shown explicitly below. But there are also ring terms 
generated in the charge sector, and the dominant one 
takes precisely the form in Eq. (1), where two Cooper 
pairs are destroyed on opposite corners of an elementary 
square plaquette, hopping to the other two corners. 

Consider then a Hamiltonian version of the Z2 gauge 
theory: 

H = H c + H a +H s , (Dl) 

H c = -t c £V; r ,(&t& r ,+h.c.) 

<rr') 



+ U^(blb r -2n) 2 , (D2) 

I" 

H a = ° Z rr> ~ h <&> > (D3) 

□ □ <rj> 

H s = -J2 °ir> [ts (Prfr' + h.C.) 

<rr') 

+ A rr , (/ rT / rU - / rU / rT + h.C.)] , (D4) 

where r and r' denote sights of a 2d square lattice. Here 
b\ creates a chargon at site r while f£ a creates a spinon 
with spin a =t,| at site r. The operator b\b T mea- 
sures the number of charge e bosons at site r, and the 
Hubbard-like U term sets the mean Cooper pair density, 
n r = &j6 r /2, to be n. With this convention, the half- 
filled electron model with charge e per site corresponds 
to n = 1/2. The constant A rr ' contains the information 
about the pairing symmetry of the spinons. The a^ r , , a^ r , 
are Pauli spin matrices that are defined on the links of 
the lattice. The full Hamiltonian is invariant under the 
Z% gauge transformation b r — > —b r , f r — > —f T at any site 
r of the lattice accompanied by letting <j* r , — > — (T^ r , on 
all the links connected to that site. This Hamiltonian 
must be supplemented with the constraint equation 

Gr = 1] v x rr ^ {lU * +blK) = 1- (D5) 

r' £r 

Here the product over a^ T , is over all links that emanate 
from site r. The operator Q T , which commutes with the 
full Hamiltonian, is the generator of the local Z2 gauge 
symmetry. Thus the constraint Q r = 1 simply expresses 
the condition that the physical states in the Hilbert space 
are those that are gauge invariant. 

Generally, the gauge field dynamics is very compli- 
cated, but it simplifies considerably when either h or K 
greatly exceeds the strength of the couplings to the mat- 
ter fields, t c ,t s and A rr /. For example, when K —> 00, 
fluctuations in a% T , are suppressed completely, and one 
can choose a gauge with the a^ r , = 1 on every link. This 
corresponds to a deconfined phase in which the spinons 
and chargons exist as bona fide particle excitations. On 
the other hand, when h — > 00 one has a^ r , ~ 1, and it 
is convenient to integrate out the gauge field perturba- 
tively. Focusing for simplicity on the Mott insulator at 
half-filling, and in the latter large h limit, with n = 1/2 
with U — ► 00, at lowest order in t s and A one recovers 
the Heisenberg spin-model, 

H 2 = J2 ]T S r -S r ,, (D6) 

<rr'> 

with J 2 = (t 2 s + A 2 )/h. Here, S r = \ft a s a pf r p, with s a 
vector of Pauli spin matrices. Note that in the h, U — * 00 
limit the gauge constraint becomes a single-occupancy 
constraint: e™^ r ^ r — 1. Upon inclusion of a small non- 
zero couplin g K, one generates the spin ring exchange 
term in Eq. ( |166| ), upon integrating out the gauge field, 
with = K(Jf+ A 4 )/(h + U) 4 . Notice that the strength 
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of the 4-spin ring exchange interaction, Jp , is propor- 
tional to the gauge theory coupling K. This suggests 
that spin models with appreciable ring exchange interac- 
tions are good candidate models to exhibit fractionalizcd 
phases. 

It is very instructive to examine the charge sector 
of the theory one generates upon integrating out the 
gauge field by similarly expanding perturbatively in small 
t c . Generally, one will generate various Cooper pair 
hopping processes. To leading order in K the Cooper 
pair plaquette ring term of Eq. ( |168| ) is generated, with 
J£=Kti/(h + U)\ 



APPENDIX E: RK MANIPULATIONS 

We begin by inserting the integral representations of 
P x / y and the explicit representation of -Pflip, inherited 
from Eqs. (148) and Eq. (143), respectively. The sums 
over {<J X y} can be explicitly performed, and we find 



L 2 components in the sum but only 2L variables). The 



(P &ip (x,y;{a})) { ^ y} = (F xy [{e^}}) {e ^ yh (El) 
where 



«™-inr(*t)° 



-«[{0*.<M] 



(E2) 

In these expressions, the transformed "flippability" func- 
tion is 



F* v [{o,<j>}]=\ n 



J x-\-u 



(E3) 



and the classical "action" is given by 



L L 

s [{QxAy}] = ~ ^ Incos(fl a - ()> y ) + imLy^(6 a - <f) a ) 

x,y=l a—1 



(E4) 



In Eq. (E4), we had added the final M 2 term to fix a 
redundancy in the description. In it's absence, the ac- 
tion s possesses a continuous translational symmetry in- 
herited from the definition of the projection operators, 

@x * @x 

+ A, (j) y — ► (j) v + A, for all x, y. This represents 
a redundancy in the constraints due to the fact that the 
total particle number of the system is given both by the 
sum of the row particle numbers and by the sum of the 
column particle numbers. The addition of the M 2 term 
has no effect on physical quantities, and has the benefit 
of lifting the unnecessary zero-mode of s. We will in fact 
take M 2 — > oo ultimately for simplicity at an appropriate 
stage of the calculation. 

No te the L prefactor in front of the second term in 
Eq. (E4), which suggests a saddle point approximation 
for large L (in the first term, note also that there are 



ds 



ds 
i ><:> 



0, are 



saddle point conditions. 

L 

J2^{6x-M = -imL- M 2 ^2{e a + <j) a ),(E5) 

y=l a 
L 

^tan(0 x -<^) = ^mL-M 2 ^(0 a + a ).(E6) 

x — l a 

These are solved by the uniform (imaginary) solution 
6 X = —c/)y = —irj/2, with 

tanhr/ = m. (E7) 

Next, expanding around this saddle, we let 6 X — ~irj/2 + 
\{ip+ x +ip- x ), 4>v = iril^ + hi^+x-^-a:)- Expanding the 
action to cubic order in ip± , as justified by the following 
analysis, one finds s = Sq + S2 + S3, with 



so 

S2 



(rn In 


1 + m 




1 — m 





(1-TO 2 ) 



S3 



+ ln(l - to 2 ) 

a 



a . b 



(E8) 



(E9) 



/ 2 
— to(1 — m 
12 v 



. (E10) 



The constant term so drops out of the observables in 
which we are interested. The quadratic action, s 2j gov- 
erns small fluctuations of ip±, treating treating cubic 
and higher order terms as perturbations. The resulting 
Gaussian averages obey Wick's theorem, governed by the 
propagators that are obtained in the usual way by invert- 
ing the quadratic form in S2- 



G ab = (lp-a1p-b)s 2 



2 / Sab 

1 - TO 2 V L 

2 / Sab 

1 — TO 2 \ L 



1 

L 2 
1 

2L 



(Ell) 



where in Eq. ( Ell ) we have tak en for simplicity the 
M 2 -> 00 limit. Eqs. ( |Ell[jE12| ) imply that the fluc- 
tuations (variance) of individual ip a ± fields are small 
(0(1 /L)), and moreover the correlations between fields 
at different a^b are even smaller (0(1/L 2 )). 

For this reason, the energy density e(m) is determined 
by the saddle-point value alone, i.e. 

(PBip(.x,y;{a})) { } = F ay [e = -<t> = -iv/2], (E13) 



leading directly to Eq. (152) 



The "single-particle" gap Ai is slightly more involved, 
since we need the energy up to terms of order 1/L. To 
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calculate it, we require Ei(m; A), which is obtained from 
the expectation value of the flippability in a state with 
A additional bosons on row and column 1. This is ob- 



tained, according to Eq. (154), by slightly modifying 
the projection operators on this row and column. This 
amounts to adding an additional source term to the ac- 
tion, s[{6>, </>}] -> s[{9, <f>}] + 2iX(0i - 4>x), but otherwise 



calculating the same expectation value as in Eq. ([El]). 
Transforming to the ip± variables as before, this addi- 
tional linear term can be removed from the action by 



shifting Va- 



2iXG al , which "completes the 



square" in the action using S2 from Eq. (E9). Depen- 
dence upon A thereby moves into F and S3: 



Ei(m;X) = -v^/F xv [6 a = -^r]+^(4) +a + i)- a -2iXG al ),(p a = ^r]+^(-il} +a 

xy » 



tf>- 



2iXGZ 



(E14) 



r 



where the prime indicates the expectation value is as de- 
fined in Eq. (E2), except that sz[tp+a, V'-a] is replaced by 



This can be eval- 
and A (the latter 



s' 3 bP+a, ip-a] = s 3 [ip +a , tp-a ~ 2iAG ol ]. 
uated by expanding F xy above in ip± 
since <C 1), and further evaluating the expectation 
value perturbatively in S3. A careful examination of these 
expansions shows that there are only two contributions 



to Ai. The first comes at 0(i[>±) from the pure saddle- 
point contribution to F xy , expanded to second order in 
A. The second comes from expanding F xy to Q( toJ a t 
A = and evaluating the expectation value in Eq. ( |E14| ) 
to first order in S3 (itself expan ded t o 0(X 2 i/j±)). Adding 
the two gives the result in Eq. (|l55|). 



1 W. W. Warren, Jr., R. E. Walstedt, G. F. Brennert, R. J. 
Cava, R. Tycko, R. F. Bell, and G. Dabbagh, Phys. Rev. 
Lett. 62, 1193 (1989). 



A. V. Puchkov, P. Fournier, D. N. Basov, T. Timusk, A. 
Kapitulnik, and N. N. Kolesnikov, Phys. Rev. Lett. 77, 
3212 (1996). 

3 A. G. Loeser, Z.-X. Shen, D. S. Dessau, D. S. Marshall, C. 
H. Park, P. Fournier, and A. Kapitulnik, Science 273, 325 
(1996). 

4 V.J. Emery and SA. Kivelson, Phys. Rev. Lett. 74, 3253 
(1995) and references therein. 

5 L. Balents, M.P.A. Fisher and C. Nayak,Phys. Rev. B60, 
1654 (1999). 

6 T. Senthil and M.P.A. Fisher, Phys. Rev. B 62, 7850 

(2000) . 

7 N. Read and B. Chakraborty, Phys. Rev. B40, 7133 
(1989). 

8 S. Kivelson, Phys. Rev. B39, 259 (1989). 

9 X.G. Wen, Phys. Rev. B44, 2664 (1991). 

10 T. Senthil and M.P.A. Fisher, Phys. Rev. Lett. 86, 292 

(2001) . 

11 T. Senthil and M.P.A. Fisher, Phys. Rev. B63, 134521 
(2001). 

12 J.C. Wynn, D.A. Bonn, B.W. Gardner, Y.J. Lin, R.X. 
Liang, W.N. Hardy, J.R. Kirtley, K.A. Moler, Phys. Rev. 
Lett. 87, 197002 (2001). 

13 J.R. Kirtley, C.C. Tsuei, H Raffy, Z.Z. Li, V G. Kogan, 
J.R. Clem and K.A. Moler, |sond-mat/0103474 . 

14 "Spin waves and electronic interactions in La2Cu04", 
R. Coldea, S.M. Hayden, G. Aeppli, T.G. Perring, CD 
Frost, T.E. M ason, S.-W. Cheong, and Z. Fisk, |cond- 
mat/0006384 



G. Misguich. B. Bernu and C. Waldtmann. Phys. Rev. 



guich, Phys. Rev. B62, 6372 (2000). 
16 "Fractionalization in an Easy-axis Kagome Antiferr omag- 



18 



net", L. Bale nts, S.M. Girvin and M.P.A. Fisher, |cond- 
mat/01 10005] . 

r ' M.P.A. Fisher, P.B. Weichmn, G. Grinstein and D.S. 
Fisher, Phys. Rev. B 40, 546 (1989). 

S. Sachdev, "Quantum Phase Transitions", Cambridge 
Univ. Press. (Cambridge, 2001) and references therein. 
M. Greiner, O. Mandel, T. Esslinger, T. Hansch and I. 
Bloch, Nature, 415, 39 (2002). 
M.P.A. Fisher, Phys. Rev. Lett. 65, 923 (1990). 
"Bond and Neel order and fractionalization in ground 
states of easy-plane antif erromagnets in two dimensions" 
K. Park and S. Sachdev, 



ond-mat/0112003 



B60, 1064 (1999); W. LiMing, P. Sindzingre and G. Mis- 



S. Tomonaga, Prog. Theor. Phys. (Kyoto) 5, 544 (1950); 
J.M. Luttinger, J. Math. Phys.N.Y. 4, 1154 (1963). 
F.D.M. Haldane, J. Phys. C14, 2585 (1981); Phys. Rev. 
Lett. 47, 1840 (1981). 

M.P.A. Fisher and L.I. Glazman, in "Mesoscopic Electron 
Transport", edited by L. Kowenhoven, G. Schoen and L. 
Sohn, NATO ASI Series E, Kluwer Ac. Publ., Dordrecht 
(1997) and references therein 

D.J. Thouless, Proc. Phys. Soc. Lond. 86, 893 (1965). 
J.M. Delrieu and J.H. Hetherington, Rev. Mod. Phys. 55, 
1 (1983). 

The duality transformation on 2d bosons that we introduce 
and exploit in this paper is explicitly tailored for bosonic 
ring models, and is quite distinct from the standard boson- 
vortex duality in 2+1-dimensions, as discussed for example 
in M.P.A. Fisher and D.H. Lee, Phys. Rev. B39, 2756 
(1989). 

AT. Dorsey, Phys. Rev. B43, 7575 (1991) 
D. Carpenti er. C. Peca, 
mat/010319^ unpublished. 



and L. Balents, cond 



28 



30 R. Shankar, Rev. Mod. Phys. 66, 129 (1994). 

31 L. Balents, unpublished (2002). 

32 The crossover to the superfluid state is slightly subtle if 
for some reason the most relevant tunneling operator has 
m > 1. In this case, these strongest tunneling processes 
conserve the boson number on m subsets of the rows and 
m subsets of the columns, leading to a state which should 
be thought of as a collection of m 2 interpenetrating and 
coexisting superfluids (i.e. with m 2 second sound modes). 
However, this state is in turn unstable to arbitrarily weak 
tunneling tn, which generates a Josephson coupling be- 
tween these superfluids. This leads to a single (fully mixed) 
second sound mode at low energies. 

33 L. Balents and M.P.A. Fisher, unpublished (2002). 

34 D.S. Rokhsar and S. Kivelson, Phys. Rev. Lett. 61, 2376 
(1988). 

35 P.W. Anderson, Science, 235, 1196 (1987). 

36 S. Kivelson, D.S. Rokhsar, and J. Sethna, Phys. Rev. B35, 
8865 (1987). 

37 PA. Lee and N. Nagaosa, Phys. Rev. B45, 966 (1992). 

38 T. Senthil and M.P.A. Fisher, page 121 in "More is Differ- 
ent: Fifty years of Condensed Matter Physics", edited by 



N. Ong and R. Bhatt (Princeton University Press, 2001). 
G. Baskaran, Z. Zou, and P.W. Anderson, Solid State Com- 
mun. 63, 873 (1987). 

G. Baskaran and P.W. Anderson, Phys. Rev. B37, 580 
(1988). 

I. Affleck and J.B. Marston, Phys. Rev. B37, 3774 (1988); 
L. Ioffe and A. Larkin, Phys. Rev. B 39, 8988 (1989); 
PA. Lee, N. Nagaosa, T.-K. Ng, and X.G. Wen, Phys. 
Rev. B57, 6003 (1998). 

See, "Breakdown of Fermi-liquid theory in a cuprate su- 
perconductor" , R.W. Hill C. Proust, L. T aillefer, P. 
Fournier and R.L. Greene, |cond-mat /0 1 12269 , and refer- 
ences therein. 

"Wigner glass, spin liquids and the metal-insulator tran- 
sition", S. Chakravarty, S. Kivelson, C. Nayak and K. 
Voelker, Philos. Mag. B 79, 859 (1999). 
"Exchange frequencies in the 2D Wigner crystal", B. 
Bernu, L. Candido and D.M. Ceperley, Phys. Rev. Lett. 
86, 870 (2001). 

M. Hermele, unpublished (2002). 



